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Abstract 

We consider a class of Hamiltonian nonlinear wave equations governing a 
field defined on a spatially discrete one dimensional lattice, with discreteness 
parameter, d = h , where h > is the lattice spacing. The specific cases we 
consider in detail are the discrete sine-Gordon (SG) and discrete (j) models. For 
finite d and in the continuum limit (d — ► oo) these equations have static kink-like 
(heteroclinic) states which are stable. In contrast to the continuum case, due to 
the breaking of Lorentz invariance, discrete kinks cannot be "Lorentz boosted" 
to obtain traveling discrete kinks. Peyrard and Kruskal pioneered the study of 
how a kink, initially propagating in the lattice dynamically adjusts in the absence 
of an available family of traveling kinks. We study in detail the final stages of 
the discrete kink's evolution during which it is pinned to a specified lattice site 
(equilibrium position in the Peierls-Nabarro barrier). We find: 

(i) for d sufficiently large (sufficiently small lattice spacing), the state of the 
system approaches an asymptotically stable ground state static kink (centered 
between lattice sites). 

(ii) for d sufficiently small d < d* the static kink bifurcates to one or more time 
periodic states. For the discrete (J) 4 we have: wobbling kinks which have the same 
spatial symmetry as the static kink as well as "g-wobblers" and "e- wobblers" , 
which have different spatial symmetry. In the discrete sine-Gordon case, the "e- 
wobbler" has the spatial symmetry of the kink whereas the "g-wobbler" has the 
opposite one. These time-periodic states may be regarded as a class of discrete 
breather / topological defect states; they are spatially localized and time periodic 
oscillations mounted on a static kink background. 

The large time limit of solutions with initial data near a kink is marked 
by damped oscillation about one of these two types of asymptotic states. In 
case (i) we compute the characteristics of the damped oscillation (frequency 
and d- dependent rate of decay). In case (ii) we prove the existence of, and 
give analytical and numerical evidence for the asymptotic stability of wobbling 
solutions. 
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The mechanism for decay is the radiation of excess energy, stored in internal 
modes, away from the kink core to infinity. This process is studied in detail 
using general techniques of scattering theory and normal forms. In particular, 
we derive a dispersive normal form, from which one can anticipate the character 
of the dynamics. The methods we use are very general and are appropriate for 
the study of dynamical systems which may be viewed as a system of discrete 
oscillators (e.g. kink together with its internal modes) coupled to a field (e.g. 
dispersive radiation or phonons). The approach is based on and extends an 
approach of one of the authors (MIW) and A. Soffer in previous work. Changes in 
the character of the dynamics, as d varies, are manifested in topological changes 
in the phase portrait of the normal form. These changes are due to changes in 
the types of resonances which occur among the discrete internal modes and the 
continuum radiation modes, as d varies. 

Though derived from a time-reversible dynamical system, this normal form 
has a dissipative character. The dissipation is of an internal nature, and corre- 
sponds to the transfer of energy from the discrete to continuum radiation modes. 
The coefficients which characterize the time scale of damping (or lifetime of the 
internal mode oscillations) are a nonlinear analogue of "Fermi's golden rule", 
which arises in the theory of spontaneous emission in quantum physics. 



1 Introduction 

Coherent structures, e.g. kinks, solitary waves, vortices, play a central role, as carriers 
of energy, in many physical systems. An understanding of their dynamical proper- 
ties, e.g. stability, instability, metastability, is an important problem. While for many 
years Hamiltonian partial differential equations and their coherent structures, defined 
on a spatial continuum, have received a great deal of attention, there has been in- 
creasing interest in spatially discrete systems. Two important reasons are that (a) 
certain phenomenoma are intrinsically associated with discreteness and (b) numerical 
approximations of continuum systems involve the introduction of discreteness, which 
may lead to spurious numerical phenomena which require recognition. Some examples 
of the use of discrete systems in the modeling of physical phenomena are: the problem 
of dislocations propagating on a lattice (for which the Frenkel - Kontorova or discrete 
sine-Gordon model was originally proposed) [^2], [I], |23|], arrays of coupled Josephson 
junctions, |55], the problem of the local denaturation of the DNA double helix 



I], [T7|, [19], [18] and coupled optical waveguide arrays |21], pE5] . 
A result of the many investigations of discrete systems over the last ten to fifteen 
years has been a recognition, mainly through numerical experiments and heuristic 
arguments, of the often sharp contrast in behavior between the dynamics of discrete 
systems and their continuum analogues. Some of these contrasts are easy to anticipate. 
For example, continuum systems modeling phenomena in a homogeneous environment 
are translation invariant in space, and may have further symmetry, e.g. Galilean or 
Lorentz invariance that enable one to construct traveling solutions from static solutions. 
The analogous discrete system is expected to lose these symmetries and therefore the 
existence of traveling wave solutions is now brought into question. A natural question 
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concerns the propagation of energy in the lattice; if one initializes the system with 
data which, for the continuum model would result in a coherent structure propagating 
through the continuum, what is the corresponding behavior for the energy distribution 
on intermediate and long times scales on the lattice? Does the system "find" a coherent 
structure to carry the energy? Does the energy get trapped or pinned? In this paper 
we seek to obtain insight into these questions for a class of discrete nonlinear wave 
equations. The special cases we consider in detail are the discrete sine-Gordon (SG) 
and discrete 4 equations. The methods however are rather general and apply to 
systems which can be viewed as the interaction between a finite dimensional system of 
"oscillators" with an infinite-dimensional system governing a continuous spectrum of 
waves; see the further discussion below and in section 7. 

The basic characteristics of the dynamics of coherent structures in discrete systems 
were systematically explored in a pioneering paper by Peyrard and Kruskal [48|; see 



also the contemporaneous papers [E7], E7L 153]. Of particular relevance to our work 



are the more recent articles of Boesch, Willis and coworkers |]59| , [54], || [5[ [11], |T(J]. In 



this work the discrete sine-Gordon equation is solved with kink-like initial data on the 
lattice. Observed is a rapid initial velocity adjustment of the kink, a quasi-steady state 
phase, resonances of the kink oscillations with phonons (continuous spectral modes) 
and the emission of radiation, which results in the kink's deceleration and eventual 
pinning in the Peierls-Nabarro (PN) potential. The final stage involves the relaxation 
to an asymptotic state, which is either a static or time-periodically "dressed" kink. 
In a rough sense, this is a result of the absence of a smooth family of traveling kink 
solutions due to broken Lorentz invariance. Insight into this process can be gleaned by 
studying the linearized spectrum about the kink. 

An important feature of the discrete systems that makes them very different from 
their continuum analogues is the presence of additional neutral oscillatory modes in 
the spectrum of the linearization about the coherent structures. The sources of these 



internal modes are principally of two types; see for example [0, |3^, |30], (1) The 
continuum system is translation invariant, a symmetry which leads to the continuum 
linearization having zero modes, eigenmodes and generalized eigenmodes corresponding 
to zero frequency. Viewed as a perturbation of the continuum system, the discrete, but 
"nearly" continuum problem is expected to have nearby modes to which these zero 
modes have been perturbed. If the coherent state is stable, these modes should be 
neutrally stable, i.e. that is, they must correspond to purely imaginary eigenvalues. 
(2) It is possible that discrete neutral modes may emerge from the continuous (phonon) 
spectrum. 

Although some of these phenomena have been identified in the early numerical 
investigations [16), there has not been a systematic dynamical systems study relating 
particular resonances to the rate with which the kink is trapped by a "valley" in the 
PN potential or, once inside such a valley, the rate with which the kink relaxes to its 
asymptotic state. The main results in this direction date from the work of Ishimori 
and Munakata PTfl , in which the McLaughlin-Scott direct perturbation scheme [II 



is 



used, valid in only in the nearly continuum regime, and the work of Boesch, Willis and 
El-Batanouny ||, which is based on numerical simulations and heuristic arguments. 
In this paper, we introduce a systematic approach to the study of these phenomena. 
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The work is based on the recent work of one of us (MIW) with A. Soffer on a time 



dependent theory of metastable states in the context of: quantum resonances [pi 



ionization type problems (parametrically excited Hamiltonians) and resonance and 



radiation damping of bound states in nonlinear wave equations |53| ; see also [33], ^2 
A fruitful point of view adopted in these works and in the present work is that the 
dynamics can be understood as the interaction between a finite dimensional dynamical 
system, governing the bound state (internal modes plus kink) part of the solution and 
an infinite dimensional dynamical system and governing radiative behavior. Using 
the tools of scattering theory and the idea of normal forms, we derive a dispersive 
normal form, which is a closed (up to controllable error terms) finite dimensional system 
governing the internal mode components of the perturbation. From this normal form 
many aspects of the the large time asymptotic state are deducible. 

The discrete systems we study {e.g. discrete sine-Gordon and discrete 4 ) depend 
on a discreteness parameter, d; as d increases the continuum limit is approached. The 
character of the normal form (topological character of its phase portrait) changes with 
d because the kinds of resonances which occur among discrete mode oscillations and 
continuum radiation (classifiable in terms of integer linear combinations of internal 
mode frequencies) change with d. 

Our analytical and numerical studies lead us to the following picture concerning 
the dynamics in a neighborhood of the kink: 

(1) The ground state kink, K gs , is always Lyapunov stable^ 

(2) If d, the discretization parameter, is sufficiently large (approaching the spatial con- 
tinuum case), the ground state kink is an attractor, i.e. is asymptotically stable. In 
this case the kink is approached at different algebraic rates depending on the range of 
d values. We infer this from the normal form analysis of section 4- 

(3) For d sufficiently small, our discrete nonlinear wave models have one or more 
branches of finite energy time-periodic solutions which bifurcate from K gs ; see Theorem 
5.1 and Corollary 5.1. For the specific cases of discrete SG and discrete <p A our results 
imply, for various regimes of d, that there exist wobbling kinks, W . These are time 
periodic solutions with the same spatial symmetry as the kink and have been previously 
observed in numerical simulations. Additionally, our results imply in certain regimes 
of d, the existence of time periodic solutions (such as gW , eW in the <p A and gW in 
the SG), u{t), with the property that u{t) — K gs is, to leading order in the direction of 
an even internal mode. Tables 3 and 1 in section 4 indicate the regimes in which these 
various solutions denoted occur. The wobbling solutions (W, gW, eW ) may be viewed 
as a class of discrete breather / topological defect states. They are spatially localized 



1 If the initial data is in a small neighborhood of K gs , then the solution remains in a small 
neighborhood of K gs for all time. This notion of stability does not however imply convergence to K gs . 

2 By asymptotic stability we mean that a small perturbation of the kink gives rise to a solution 
which converges as t — > ±oo, in some physically relevant norm, to a kink. Asymptotic stability is 
a notion of stability commonly associated with dissipative dynamical systems and is not commonly 
associated with general energy conserving systems. In this work we are studying infinite dimensional 
Hamiltonian systems on infinite spatial domains. These systems have the possibility of radiation of 
energy to infinity, while keeping the total energy of the system preserved. Thus dissipative behavior is 
realized through dispersion and eventual radiation of energy out of any compact set; see, for example, 
the short overview in fl57fl and references therein. 
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and time-periodic oscillations mounted on a static kink background. The usual discrete 
breathers || ; [J<|/ are mounted on a zero background. This is analogous to the situation 
with solitons. Dark solitons of the defocusing nonlinear Schrddinger (NLS) equation 
are mounted on a non-zero continuous wave background, while standard solitons of 
focusing NLS sit on a zero background. 

(4) Numerical simulations and the normal form analysis indicate that in various d 
regimes these time periodic states are local attractors for the dynamics. 

(5) The analysis of section 5 implies some nonexistence results concerning solutions 
of discrete nonlinear wave equations in a neighborhood of the kink. We have that in 
certain regimes of the discreteness parameter, d, no nontrivial time-periodic solution 
exists, and that time quasiperiodic solutions do not exist in any regime of d. On the 
other hand, the normal form analysis and numerical simulations indicate, in some 
regimes of d, that periodic or quasiperiodic oscillations can be very long lived; see 
section 6, and in particular, the discussion of Regime VI in section 6.1. Thus, we may 
think of the system as possessing metastable periodic and quasi-periodic solutions. In 
regimes of d where we show that the wobbling kink, W , is unstable on long time scales 
due to a resonance of the "shape mode" with continuum radiation modes, we have the 
analogue of the the wobbling solution of continuum <f> A , proved to be stable on large but 
finite time scales by Segur ]56|/ . On an infinite time scale this wobbling solution behaves 
as those of the discrete system for large enough d; eventually the oscillations damp at 
an algebraic rate leaving a kink in the limit flUQ , \4QJ . For the continuum system, the 
limiting kink may be traveling, while for the discrete system it is a static ground state 
(centered between lattice sites) kink. 

In (2) the dispersive normal form has a dissipative character; the effect of coupling 
the internal mode oscillations to radiation is modeled by an appropriate nonlinear fric- 
tion; see section 4 and the discussion in the introduction to [[53] on another related 
model. The damping or friction coefficients are given by formulae which can be un- 
derstood as a nonlinear generalization of Fermi 's golden rule, arising in the context of 
the theory of spontaneous emission in atomic physics. Finally, it is worth noting that 
we obtain a dissipative and therefore apparently time-irreversible normal form from a 
system which is conservative and time-reversible. There is no contradiction because the 
dissipation is of an internal nature; it signifies the transfer of energy from the discrete 
internal mode oscillations to the continuum dispersive waves which propagate to infin- 
ity. That dissipative dynamical systems emerge from conservative systems which are 
a coupling of a low dimensional dynamics to infinite dimensional dynamics ("masses 
and springs coupled to strings") has been observed in many contexts; see, for example, 
3~5| , ||, |37| and the discussion in [52], [53| . 



The paper is organized as follows: 

• Section 2 begins with a general discussion of discrete nonlinear wave equations 
which support kink-like (heteroclinic) structures and then specializes to a discus- 
sion of the discrete sine-Gordon (SG) and discrete </> 4 systems. 

• In section 3, we present the decomposition of solutions into discrete internal mode 
components and radiative components, enabling us to view the dynamics near a 
kink as the interaction of finite and infinite dimensional Hamiltonian systems. 
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In section 4 the dispersive normal forms are derived and discussed for both dis- 
crete SG and discrete 4 . 



In section 5, we prove that if f2 is an internal mode frequency and no multiple 
of it lies in the continuous spectrum of the kink (phonon band) , then there is a 
family of finite energy time periodic solutions which bifurcate from the kink in 
the "direction" of the corresponding internal mode. The proof is based on the 
Poincare continuation method and is an application of the implicit function the- 
orem in an appropriate Banach space. This approach was used also by MacKay 



and Aubry []38|], who constructed discrete breathers of nonlinear wave equations 



in the anti-integrable limit. 

In section 6 we combine the normal form analysis of section 4 and existence 
theory for periodic solutions of section 5 with observations based on numerical 
simulations to obtain a detailed picture of the dynamics in a neighborhood of the 
ground state kink. 

Finally, in section 7, we summarize our results and give directions of interest for 
future research. 



1.1 Notation 

• 7L denotes the set of all integers and IR denotes the set of real numbers. 

• For u = {wi}i e z, 8\u denotes the discrete Laplacian of u defined by: 

{5 2 h u)i = hr 2 (u i+ i - 2u,i + Uj_i) . (1.1) 
In the case h = 1 we shall use the simplified notation 5 2 = 8f. 

• The inner product of vectors u and v in 1 2 {7L) is given by 

(u,v) = ^UjVj. (1.2) 

j 

• The t subscript denotes time partial derivatives, e.g. u t (t) = d t u(t), whereas i 
denotes the lattice site numbering. 

• 1 2 {7L) is the Hilbert space of sequences {Mj}i e z which are square summable. 

• H S (I) denotes the Hilbert space of functions /, defined on / C IR such that / and 
all its derivatives of order < s are square-integrable. For the case of 2ir periodic 
functions, / = S^, an equivalent norm on H s is given by: 

\\f\\h = 'L( 1 + \ n \ 2 Y\M 2 > (i-^) 

where f n denotes the n th Fourier coefficient of /: 

/„ = (27r)-3 r e~ 2mnt f(t) dt. (1.4) 



o 
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• For example, if 2 (IR; 1 2 {7L)) denotes the space of functions f(t,i) which are H 2 , 
as functions of t, with values in the space of 1 2 {7L) functions. 

• cr(B) and a cont (B) denote respectively the spectrum and the continuous (phonori) 
spectrum of the operator B. 

• < d denotes the discretization parameter, d large is the spatial continuum 
regime, while for d small the effects of discreteness are strong. 

2 Discrete and continuum wave equations 
2.1 General Background 

Two model nonlinear wave equations on lattices that have played a central role in the 
theory of nonlinear waves are the discrete sine- Gordon equation (SG) and the discrete 
4 model (0 4 ). These equations may be viewed as governing the dynamics of a chain 
of unit mass particles which, in equilibrium, are equally spaced, a unit distance apart. 
The particles are then subjected to a conservative force derived from a potential. In 
terms of the displacement of the i th particle from equilibrium the equations of motion 
areQ: 

Ui, u = (S 2 u)i - d- 2 V'( Ui ). (2.1) 
The case of SG and 4> 4 correspond, respectively, the choices of potential: 

V{u) = 1-cosu, (SG) 
V(u) = \{l-u 2 ) 2 . (0 4 ) 

More generally, we assume V has three continuous derivatives and satisfies certain 
constraints appearing below, related to the existence of kink solutions or homoclinic 
orbits. 

The parameter d is a fixed constant with, d~ 2 , having the interpretation of the ratio 
of on- site potential energy to elastic coupling energy. 

Relation between discrete and continuum systems: The parameter d can also be 
seen to play the role of the reciprocal of the lattice spacing in passing between contin- 
uum and discrete models. Consider the continuum model governing the displacement 
v (T , xj '. 

vtt = v xx -V\v). (2.2) 

Introducing the lattice spacing parameter, h, and replacing u xx (x) by (S 2 u)i we obtain 
the discrete nonlinear wave equation governing Vi = v(T,i; h): 

v l ,TT = (S 2 h v) l -V'(v l ). (2.3) 

The form of the discrete nonlinear wave equation we use is obtained by setting t = hr x T 
and defining the time-scaled displacement u(t, i; h) = v(T, v, h). Then, u(t, i; h) satisfies 
(p|) with d = h-\ 

3 We use the notational conventions introduced by Peyrard and Kruskal Eq|. 
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Taking the formal limit d j oo or equivalently h { we have that 



U(T, i) = lim u(dT, z; cT 1 ) = lim w(T, i; cT 1 ) (2.4) 



a!— >oo 



satisfies the associated continuum nonlinear wave equation ( |2.2|) . 
Hamiltonian structure: The discrete and continuum equations we consider are in- 
finite dimensional Hamiltonian systems with Hamiltonian energy functionalsQ: 

W = E + 5 (^+1 - + (2.5) 

i 

in the discrete case and 

W = / tt(«t + M D + ^0) ( 2 - 6 ) 

iR 2 

Static kink solutions: Consider time-independent or static solutions of these dy- 
namical systems. In each case, such solutions satisfy the equation obtained by setting 
time derivatives equal to zero. Thus, in the discrete case we have: 

(S'K), = V'{K t ), (2.7) 

and in the continuum case 

K"(x) = V'{K{x)) (2.8) 

Spatially uniform solutions occur at the critical points of the potential, V. Of particular 
interest are the static kink solutions. These are static solutions which are heteroclinic. 
That is, they are static solutions which can be viewed as connections, as i — > ±oo 
(respectively x — > ±oo) in the phase space of two distinct "unstable" equilibria, values 
Kf for which 

V(K+) = V(K~), V'(Kt) = 

V" = V"(Kt) = V"(K;)>0. (2.9) 

Kink solutions can also be constructed by variational methods. Consider the Hamil- 
tonian energy functional, Tt restricted to t-independent functions: 

M M ] = ^ f 2^ Ui+1 ~ Ui ^ 2 ~^ d~ 2 V( u i) ^ ? discrete case (2-10) 
h[u] = [ -u 2 + V(u) dx continuum case. (2-11) 

JR 2 

Then, the Euler Lagrange equations associated with h[u] is the equation for the kink, 
( j2.7p . That is, if K gs denotes a minimizer of h[u] then since for all ip G / 2 (/Z), 

— h[K gs + T^}\ T=0 =0. (2.12) 



4 There is freedom in the choice of potential V; we choose V so that the Hamiltonian energy of 
the static kink is finite. 
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Figure 1: Ground state static kink for the very strongly discrete sine-Gordon model 
(-o-) d = 0.6, and very close to the continuum limit d = 10. 



This implies that for all ip G l 2 (TL) 

E [-(S 2 K 9 s)i + d- 2 V'{K, 



0. 



(2.13) 



This is equivalent to (|2.7|) . A solution constructed by minimization of h[u] is called a 
ground state kinkn 

Invariance and broken invariance: 

An important structural difference between the discrete and continuum cases is 
that the continuum case has greater symmetry. In particular, ( |2.2j ) has the property 
of translation and Lorentz invariance, while ( |2.1|) does not. A consequence of this is 
that static solutions of (I272T) can be translated: 



K(x) i — > K(x — xq) 



(2-14) 

5 A proof that the minimum of the functional ( |2.10|) is attained can be given using the following 
strategy, used in a similar problem |p^| : For positive integers, N, define h^iu] to be the truncated 
Hamiltonian, where the summation is taken over —N < n < N. Consider h^, restricted to vectors 
satisfying u±n — K^. For any admissible vector, u = {ui}^ <N , it is possible, by rearrangement, to 
replace it with another, v, which is monotonically increasing and for which h[v] < h[u]. Therefore the 
minimizer of /ijvH; Kgs , is monotonically increasing from K~ to . A priori estimates, derived 



from the Euler-Lagrange equation (2.7), and the boundary conditions enable one to pass to the limit 



and obtain a minimizer of h[-\, K gs = {K gs ^}i e z, a- monotonically increasing static kink. 
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to give a recentered kink and Lorentz boosted: 



K(x) i-> K((x - vt)/y/l - v 2 ), \v\ < 1, (2.15) 

to give a traveling solution of velocity v. These symmetries are absent in the discrete 
models. Discrete models do however have a discrete translation symmetry but in the 
models considered this does not give rise to discrete traveling waves. 
Dynamic stability of kinks: The issue of dynamic stability is a subtle one. First, 
there is the question of what is the object we expect to be stable. Also, since our 
systems are infinite dimensional and not all norms are equivalent, different choices of 
norms will measure different phenomena. Here we briefly discuss three related notions 
of stability we have in mind: (a) orbital Lyapunov stability, (b) asymptotic stability 
and (c) linear spectral stability. 

(a) Orbital Lyapunov stability means stability of the shape of the kink; if at t = the 
data is nearly shaped like a kink then it remains shaped like a kink for allt ^ 0. The 
kink-like part of the solution will typically move so the solution at different times is 
near some time-dependent spatial translate of the static kink (element of the symmetry 
group orbit of the static kink). 

In the continuum case, orbital Lyapunov stability in the space H 1 (1R) is a conse- 
quence of characterization of the kink as a local minimizer of 7i. A detailed proof is 



presented in [J^]. A simple proof of Lyapunov stability of ground state discrete kinks 
can be given based on the same ideas. 

(b) Orbital asymptotic stability means that if the initial data is nearly a kink then the 
solution converges to kink (possibly traveling, in the continuum case or pinned in the 
discrete case) as t — > oo. 

For both notions (a) and (b) we see that the stable object is the family of solitary 
wave solutions. That is, in the case where there is a multiparameter family of solutions 
(kinks, solitary waves ...), it is the collection of all such that is the stable object. 
Concerning the terminology orbital stability, this family of solutions is related to the 
set of functions generated by applying elements of the equation's symmetry group to 
the solution (here, static kink), and thus we can think of the group orbit as stable. 

(c) Spectral stability: 

As is well known, dynamic stability of the particular solution is related to the 
spectrum of the linearization of the dynamical system about this solution. Lineariza- 
tion about the kink gives an evolution equation for infinitesimal perturbations of the 

6 The proof is based on the following idea. Since K gs is an energy minimizer, by fl2.24| ), the second 
variation at the kink, B 2 is non-negative. In fact, it is strictly positive; that is, for all -0 € ^ 2 (^)> 
(tp,B 2 ip) > uig(ip,i/j), where lo 2 — uig(d) > is the discrete eigenvalue to which the zero mode of 
the continuum equation perturbs due to discretization of space. Suppose we have initial data near a 
kink: u(0) = K gs + Vq, d t u(0) = v\ and let u{t) = K gs + v(t) be the resulting solution. Assume that 
11^0) fill; 2 ~ £ : where e is sufficiently small. Then, since K gs is a critical point of h[-\ we have: 

E~h[K gs +v{t)]\ t = -h[K gs ] = h[K gs +v(t)]-h[K gs ] = ^{v(t),B 2 v(t))+0(\\v(t)\\%) 

>^ll«(*)llp-c||«(*)||? a . 
implying that ||w(t)||p ~ 0{e) for all t ^ 0, and K gs is stable. 
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following form: 

For the discrete case: 



(d 2 + B 2 )p = 0. (2.16) 



B 2 Vl = -(5 2 p) l + d- 2 V"{K t ) Pi) (2.17) 
and for the continuum case: 

B 2 p = -p xx + V"(K(x))p. (2.18) 

At a minimum we expect that stability can hold only if no solutions of the linearized 
evolution equation, corresponding to finite energy initial conditions, can grow as \t\ 
increases. In carrying out linear stability analysis we seek solutions of the form pi = 
e xt Pi (respectively, p = e xt P(x)) and obtain linear eigenvalue problems of the form: 



P 2 + A 2 )P = 0. (2.19) 



Explicitly, we have: 



- A 2 P = -(5 2 P)i + d~ 2 V"(Ki)P u discrete case 
— A 2 P = —P xx + V"(K(x))P. continuum case 

We say A is in the I 2 spectrum (respectively, L 2 spectrum) of the linearization about 
the kink, or simply spectrum of the kink if the operator B 2 + A 2 does not have a bounded 
inverse on I 2 (respectively L 2 ). The spectrum will typically consist of two parts: (i) the 
point spectrum consisting of isolated eigenvalues of finite multiplicity, for which the 
corresponding solution of ( |2.19| ) is in the Hilbert space, and due to the infinite extent of 
the spatial domain (ii) continuous spectrum, whose corresponding solutions (radiation 
or phonon modes) of Q2.19D are bounded and oscillatory over the entire spatial domain. 



Note that the eigenvalue equation ( |2.19| ) has the following symmtery: if A has the 
property that ( |2.19| ) has a nontrivial solution solution then —A, A and —A also have this 
property. To avoid exponential growing solutions, we must require that the linearized 
spectrum is a subset of the imaginary axis. If this holds we say the solution is spectrally 
stable. 

Spectrum of the kink: 

We are interested in the set of A's for which the eigenvalue equation 

(P 2 + A 2 )P = (2.20) 

has a nontrivial solution. We begin by discussing the spectrum of B 2 and then take 
square roots to a obtain the spectrum of the the linearization about the kink. 

We now make two general remarks about the spectrum of the kink solution, one 
concerning the continuous spectrum and one concerning the point spectrum. 
Continuous spectrum: The continuous spectrum is determined by the solutions of the 
constant coefficient equation obtained from ( |2.19| ) by evaluating its coefficients at spa- 
tial infinity. In the discrete case we obtain: 

- A 2 P = -(5 2 P) i + d~ 2 V:P u (2.21) 
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where V" = \im i ^± 00 V"(K i ). This constant coefficient equation can be solved in 
terms of exponentials which yield the character of the exact continuum eigensolutions 
of ( p,19| ). We let A = iu. We find bounded oscillatory solutions of the form: P n = 
exp(v / — lkn), n e 7L, k real and arbitrary, where satisfies the dispersion relation: 

uo 2 = 4 sin 2 (A;/2) + d~ 2 V". (2.22) 

It follows that the continuous spectrum of B 2 is the positive interval from d~ 2 V" to 
d~ 2 V" + 4. Therefore, the continuous spectrum of the linearization about the kink 
solution consists of two intervals on the imaginary axis: ±z[d _1 y/V£, \J^ + d~ 2 V"]. 
Point spectrum: An important conclusion, concerning the point spectrum of kink for 
the continuum equation, can be made as a consequence of the equation's symmetry 
group. If K(x) is a kink then translation invariance implies that K'(x) is a zero mode, 
a solution of the linear eigenvalue problem with eigenvalue zero, B 2 K = 0. This 
zero mode is often called the Goldstone mode. Since the eigenfunction, K' does not 
change sign, zero is the ground state energy (lowest point in the point spectrum) of 
the Sturm-Liouville operator — d 2 + V"(K(x)). This eigenvalue is of multiplicity one. 

If one views the discrete model as a perturbation of the continuum model, due 
to the absence of corresponding symmetries in the discrete problem one expects the 
eigenvalue at zero to move, as d decreases from infinity (h, the lattice spacing, increases 
from zero), to the right of zero or to the left of zero. In terms of the spectrum of the 
kink (plus or minus the square root of the spectrum of B 2 ), this means that in the 
discrete case the spectrum of the kink consists of either (a) a purely imaginary pair, 
±iu s (stable case) or (b) a pair of real eigenvalues, symmetrically situated about the 
origin (unstable case). We shall see both cases in our study of specific discrete models 
and we shall refer to these eigenvalues and corresponding modes loosely as Goldstone 
modes as well. 

As we shall see in the specific models considered, other eigenvalues may appear 
in the gap between the upper and lower branches of the continuous spectrum. Such 
modes of the linearization which lie in this gap are called internal modes, and their 
number and location can change with d. For SG we shall find that there are one or two 
pairs of internal modes, while for the 4 model there are two or three pairs of internal 
modes. 

Finally, we remark that in the case of a ground state kink, one obtained by mini- 
mization of h[u], the spectrum of B 2 is non-negative and the kink is spectrally stable. 
To see this, let K gs denote a minimizer of h[u], which as indicated above, is a critical 
point of h[-]. Furthermore, for all if) G 1 2 {7L), 

d 2 

— h[K gs + rip] | r=0 > 0. (2.23) 

A calculation yields that this is equivalent to: 

~{if),B 2 ip) > 0, for all ip e 1 2 {7L). (2.24) 

Therefore, d 2 h[K gs ] = \B 2 , the second variation of h[-] at the ground state kink, is a 
nonnegative self-adjoint operator. It follows that the spectrum of the kink, ±ia(B), 
lies on the imaginary axis. 
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2.2 The Sine-Gordon equation 

Consider the discrete SG equation: 

Ui,tt = {u i+ i + Ui_ x - 2ui) - sin Uj. (2.25) 

A static discrete 2n kink is a time-independent solution, { A",},,- z of ( |2.25|) which sat- 
isfies the boundary conditions at infinity: 

Ki — > 0, as i — > — oo, and Aj — > 27r, as z — >■ oo (2.26) 

It is well-known that there are two static kink solutions, a high energy one centered 
on a lattice site and a low energy one centered between two consecutive lattice sites 
27] , [15], ^9], 0, ||, ||, [TT], 0, [34], [30], [31]]. The low energy kink corresponds to the 



minimizer of the static Hamiltonian energy, h[u], displayed in ( |2.10|) . As d increases 
(the continuum limit) the energy difference between the two static kink solutions, the 
so-called Peierls-Nabarro (PN) barrier, tends to zero and the scaled limit u(i; d~ l ) = Aj, 
with i/d = x fixed, converges to the continuum SG static kink solution 

Ksg(x) = 4tan~\e x ). (2.27) 

Linear stability analysis about the high energy and low energy kinks was carried out 
in ||. From the general discussion of section 2.1, the dispersion relation defining the 
continuous (phonon) spectrum is uj 2 = d~ 2 + 4sin 2 (/c/2). Therefore, B 2 has a band of 
continuous spectrum of length 4, [d~ 2 ,A + d~ 2 }, and the continuous spectrum of the 
kink consists of the two bands, ±i[d _1 , yA + dr 2 ] . 

For the case of the high energy kink, the point spectrum of B 2 contains a negative 
eigenvalue derived from the zero (Goldstone) mode associated with the continuum 
(translation invariant) case. Since the eigenvalue parameter in ( |2.19| ) is —A 2 , it follows 
that the spectrum of the high energy kink (parametrized by A) consists of two discrete 
real (Goldstone) eigenvalues, one positive and one negative, ±uj us . These occur at a 
distance of order exp(—n 2 d) ^Z7^ [|, pSfl . 

The low energy kink, minimizes the Hamiltonian. The second variation of h[u] 
at the kink, B 2 , is therefore a self-adjoint and nonnegative operator. Its spectrum is 
therefore nonnegative and so the spectrum of the kink is purely imaginary. In this case, 
the Goldstone mode associated with translation invariance of the continuum problem 
gives rise to point spectrum consisting of a complex conjugate pair of eigenvalues 
with order of magnitude exp(— ir 2 d). Furthermore, for d > d e , d e ~ 0.515 a spatially 
antisymmetric edge mode of energy uj 2 (B 2 e = uj 2 e) appears whose corresponding 
energy has a distance from the phonon band edge of order 0(d~ A ) 1 for d large |34], |3T| . 
This mode is not present in the continuum SG. 

2.3 Discrete 4 Model 

Consider the discrete 4 equation: 

Ui,tt = (u i+ i + Ui_i - 2ui) + ^{ui~ u 3 i) (2.28) 
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Spectrum of discrete SG ground state kink 
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Figure 2: Upper panel: schematic of spectrum of B 2 , a(B 2 ), consisting of two (for 
d > d e ~ 0.515) positive discrete eigenvalues (cu 2 < cu 2 ) and a finite band of continuous 

spectrum extending from = 1/d 2 to uj\ 
of spectrum of the kink, given by ±ia(B). 
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Discrete SG ground state kink and its internal modes 
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Figure 3: Three panels in order correspond to (a) discrete SG kink, (b) spatially even 
Goldstone mode (c) spatially odd edge mode (present for d > d e ~ 0.515) 
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In most ways the situation is quite similar to that discussed for the discrete SG. 
The 4 model has low and high energy kinks. The low energy kink is centered between 
lattice sites while the high energy kink is centered at a lattice site. As d increases 
the energy difference, the PN barrier, tends to zero, and finally, in the scaled limit 
u(i; d~ l ) = Ki, with i/d = x fixed, converges ot the continuum 4 kink, 



K^(x) = tanh(x/V2). (2.29) 

Most qualitative properties of the spectrum of the linearization about the 4 kink 
are analogous to those in the SG case. The single key difference can be traced to the 
spectrum associated with the continuum case. In addition to the zero (Goldstone) 
mode, the spectrum of the operators B 2 has internal mode of odd parity with eigen- 
frequency u 2 = 3 /(4c? 2 ). Therefore, the spectrum of the kink has purely imaginary 
internal modes associated with the energies uj s = ±iy/3/(2d). These modes are called 
shape modes. 

Turning now to the discrete 4 model we see that B 2 for the discrete case has the 
following properties: 

• The dispersion relation is uj 2 = 2/d 2 +4 sin 2 (/c/2), and therefore B 2 has continuous 
spectrum extending from 2/d 2 to 2/d 2 + 4. 

• B 2 has a positive internal Goldstone mode with energy uj 2 , so that B 2 g = ou 2 g. 
This mode, which is traceable to the zero mode of the continuum model, is 
spatially even and u g is exponentially small in d 



B 2 has an internal odd shape mode with energy uj 2 , so that B 2 s = uj 2 s. This 
mode, which is traceable to the internal shape mode of the continuum <p A model, 
is spatially odd. 

For d > d e (d e ~ 0.82) B 2 has an internal edge mode which is even and with 
corresponding energy u 2 , with B 2 e = u 2 e, which emerges from the continuous 



spectrum. For large d we have p4l pi 



j 2 ~ 2d' 2 - ^-d-\ (2.30) 



It follows that the spectrum of the discrete kink consists of pair of two or three 
complex conjugate pairs of internal modes at energies ±iu) g , ±iu s , and for d > d e , 

In the next section we shall see that whether or not certain integer linear combina- 
tions of the internal mode frequencies land in the continuous spectrum determines the 
large time asymptotic behavior. 
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Spectrum of discrete </> 4 ground state kink 
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Figure 4: 

d > d P r 



Upper panel: schematic of spectrum of B 2 , a(B 2 ), consisting of three (for 
' 0.82) positive discrete eigenvalues {u 2 < u 2 < u 2 ) and a finite band of 



continuous spectrum extending from u^m = 2/d 2 to uj^^ = ^A + 2/d 2 . Lower panel: 
schematic of spectrum of the kink, given by ±io~(B). 



3 Resonances and radiation damping 

That we can Lorentz boost a static kink and obtain a traveling kink on the spatial 
continuum (section 2) motivates the following question: 

What happens if we attempt to cause a kink to propagate through the lattice? 

The numerical and formal asymptotic investigation was initiated in ||48|| , [ 47 1 . For 
example, choose as initial conditions for the discrete sine Gordon system the state 
Ui(t = 0) = K(i/\/l — v 2 ), where K(x) is the continuum kink and < \v\ < 1. Disper- 
sive radiation plays an important role in the dynamics. The kink moves approximately 
as particle under the influence of the sinusoidal (Peierls-Nabarro) potential. As the 
kink propagates, its velocity alternately increases and decreases, depending on the lo- 
cation of its center of mass relative to the peaks and troughs of the potential. This 
oscillation leads to a resonance with the continuous spectrum j|8], [8| and a transfer of 
energy from the propagating kink to radiation modes. Eventually, the kink has slowed 
so that its energy no longer exceeds the PN barrier and then executes a damped oscil- 
lation about a fixed lattice site to which it is pinned. It is this latter process that we 
study here. 
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Discrete </> 4 ground state kink and its internal modes 
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Figure 5: Top left panel as in figure one but for 4 and d = 0.9. Top right shows the 
g eigenf unction, bottom left the Si and right the e,. 



The goal of this study is to, by analytic and numerical means, 

• identify what the asymptotic state of the system is in different regimes of the 
discretization parameter d, 

• find the rate of approach to this state. 



We begin in the setting of the discrete <f) A model, Q2.28Q . This case is a more 
complicated than discrete SG; discrete <j) 4 has either two or three internal modes, while 
discrete SG which has at most two. After a detailed treatment of the 4 case, we'll 
give an outline of the analogous results for SG. 

We begin with a solution which is a small perturbation about a low energy (ground 
state) kink. This corresponds to solving the initial value problem for Q2.28j ) with initial 
conditions: 

Ui(0) = K { + ev 0i , d t v,i(0) = ev u , (3.1) 
where e a small paramater. We then view the solution as 



Ui(t) = Ki +evi(t), 



(3.2) 



where Vi(t) denotes the perturbation about the kink. Equation fl2.28|) implies an evo- 
lution equation for vf. (df + B 2 ) v = ... and since the data for Vi is small, it is natural 
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to expand Vi in terms of the modes of the linearization about the kink. Recall that 
the operator B 2 is positive and self-adjoint, and has, depending on the range of d's 
considered, two or three internal modes g, s, and e (if d > d e ): 

B 2 g { = u 2 g g h B 2 Si = u 2 s s u B 2 ^ = uj 2 e 6i 
< u 2 < u 2 < u 2 e < 4 + 2c/~ 2 , 

where 

B 2 ^ = -5 2 ^ + d~ 2 {3K 2 - (3.3) 

We assume these internal modes to be normalized so that the set {g, s, e} is orthonor- 
mal in l 2 (TL). We also introduce the operator which projects onto the orthogonal 
complement of the subspace of internal modes: 

P c ip = ip - (g,ip)g - (s,ip)s - (e,i/j)e (3.4) 

P c is the projection onto the continuous spectral part of B 2 (phonon or radiation 
modes) and therefore the solutions of (d 2 + B 2 )u = with data in the range of P c are 
expected to decay dispersively as t — ► ±oo. 

We expand the perturbation evi in terms of the internal mode subspace and its 
orthogonal complement and thus 

Ui(t) = Ki + ea(t)gi + eb(t).Si + ec(t)ei + £ 2 ^(t), (3.5) 

with 

(g,r l (t)) = (s,r l (t)) = (e,r ] (t)) = 0, 
PcV = r] 

Substitution of Q into (ggp yields 

(a tt + u g 2 a)gi + (b tt + u s 2 b)si + (c u + uj e 2 c)ei + ed 2 r)i 

= -eB 2 r]-d- 2 (3eKiM 2 + e 2 Mf + Qe 2 K i r ] M i + 3e 3 M 2 r ]i + 0(e 4 )) (3.6) 

where: 

Mi^atygi + btysi + ctyei. (3.7) 

Notice that from here on for simplicity and compactness we drop the spatial indices i. 

We next project (|3.6|) onto the internal modes g, s and e, as well as on the range 
of P c . This yields the following coupled system of four equations for the three internal 
mode amplitudes and the continuous spectral part: 

d 2 a + u 2 a = --^(3e(g,KM 2 )+e 2 (g,M 3 ) + 6e 2 (g,KMr ] ) + 6e 3 (g,M 2 r ] ))+0(e 4 ) 

(3.8) 

d 2 b + u 2 b = -^(3e(s,KM 2 ) +e 2 (s,M 3 ) + 6e 2 (s,KM7]) +6e 3 (s,M 2 r])) + 0(e A ) 

(3.9) 

d 2 c + u e 2 c = -— (3e(e,KM 2 ) +e 2 (e,M 3 ) + Qe 2 (e,KM7]) + 6e 3 (e, M 2 V )) + C(e 4 )) 

(3.10) 

d 2 r] + B 2 r] = ~^Pc (eM 3 + 3KM 2 + 0(er]M) + 0(e 4 rj 3 )) ■ (3.11) 
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The initial data for a(t), b(t), c(t) is: 



o(0) 
6(0) 
c(0) 
77(0) 

For simplicity we shall assume: 



(g,v ), d t a{0) 

(S,V ), d t b(0): 

(e,v ), d t b(0) . 
P c v , d tV (0) = 



= (g,vi) 
= (s,vi) 
-- (e,«i) 

PcVl 



r/(0) = 0, <9 t r/(0) = 



(3.12) 



(3.13) 



corresponding to perturbations of the kink which only excite the internal modes. The 



general case can be treated as well; see |53| . 

The system ( |3.8| - |3.1lD may be viewed as finite dimensional Hamiltonian system gov- 



erning three oscillators with amplitudes a, b and c and natural frequencies ±w g , ±u s , ±u; e , 
coupled by nonlinearity to an infinite dimensional Hamiltonian wave equation for a 
field 7/. Systems of this type have been analyzed rigorously and the behavior of their 
solutions determined on short, intermediate and infinite time scales |52] , |53| . The 
analysis we present uses and extends the methods of these works. Roughly speaking 
we transform the system ( p.8|) -( ^TTD into an equivalent dynamical system, which is 
a perturbation of a normal form for which one obtains information on: (i) what the 
asymptotic state of the system is, and (ii) how this asymptotic state is approached. 



4 The dispersive normal form 

Our goal in this section is to derive the normal forms which give detailed information 
on the dynamics in a neighborhood of the ground state static kink of SG and 4 . In 
particular, the normal forms anticipate the existence of time-periodic solutions. These 
predictions are upheld by the rigorous results of section 5. In section 6, the normal form 
analysis is then joined with the existence theory of section 5 and numerical simulations 
to fill out the picture of what the dynamics are in a neighborhood of the kink. 



4.1 The normal form for discrete 4 

At this stage, we have in (|3.8| - |3.1ll) a formulation of the discrete 4 dynamical as a 
system governing the discrete internal modes interacting, due to nonlinearity, with a 
system governing dispersive waves (phonons). Our goal in this section is to present 
and implement a method, based on |52| , |53| , leading to a reformulation of the coupled 
discrete-continuum mode system (|3.8| - |3.1ll ) as a perturbed dispersive normal form in 



which the nature of the energy transfer among modes is made explicit. 

Ideally, one could solve the equation for the dispersive part, 7/, explicitly in terms 
of the discrete mode amplitudes: a, b, and c and then substitute the result into the 
mode amplitude equations to get a closed system of equations for a, b, and c. This then 
could be used in the equation for 7/ to determine its behavior. Due to the system's 
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nonlinearity, one can't expect to solve for the exact behavior of r\ in terms of a, b and 
c. Instead we solve perturbatively in e. Let 



r} = J2 £3r l 

j=0 



(4.1) 



Then, 

<9 t V°) + B 2 ri<® = -3d' 2 KM 2 . (4.2) 

The function rj captures the leading order radiative response, to the internal mode 
excitations. 

Note that r/^ is the solution of a forced wave equation. The forcing term in ( fOJ ) 
involves M 2 , which by ( p.7|) , contains squares and cubes of a, b and c. Also, note that 
since the nonlinear coupling terms on the right hand sides of equations (3.8- 3.10|) are 
small, a(t),b(t) and c(t) are slow modulations of the exponentials e ±tuj!)t , e st and 
e ±tuJet . The forcing on the right hand side of ( |4.2| ) will be resonant if nonlinearity, when 
acting on a, b and c, generates frequencies which lie in the spectrum of the operator B. 
Since B has a band of continuous spectrum, it can easily happen that integer linear 
combinations of the internal mode frequencies can lie in the continuous spectrum of B. 

Figure 6 in section 4 displays the locations of the frequencies u g ,u s and u e , their 
multiples and certain integer linear combinations relative to the continuous (phonon) 
spectrum. Figure 7 is the analogous plot for discrete SG. 

The key calculation is to compute the effect of such resonances which result in the 
transfer of energy from the discrete to continuum modes. 

Let 



a(t) = A(t) exp(iu g t) + A(t) exp(— iu g t), 
b(t) = B(t) exp(iu s t) + B(t) exp(— iu> s t), 
c(t) = C{t) exp(ici; e t) + C{t) exp(— iu> e t), 

We further impose the constraints: 

A t exp(iuj g t) + At exp(-iugt) 
B t exp(iuj s t) + B t exp(—iLJ s t) 
Ct exp(iu) e t) + Ct exp(— iuj e t) 



(4.3) 
(4.4) 
(4.5) 





0. 



and obtain: 



A t 



* t F 1 

2 



where Fa 



(2iujgj 

(2too e )- 1 e-^ t F 3 
Fn + Fi2 (i=l,2,3). The 3x2 matrix Fij is given by 



(4.6) 
(4.7) 
(4.8) 



-dr 



3e(g, KM 2 ) + e 2 (g, M 3 ) 6e 2 (g, KMt]) + 3e 3 (g, M 2 7]) 
3e(s, KM 2 ) + e 2 (s, M 3 ) Qe 2 (s, KMr\) + 3e 3 (s, M 2 rj) 
3e(e, KM 2 ) + e 2 (e, M 3 ) 6e 2 (e, KMrj) + 3e 3 (e, M\) 



(4.9) 
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We first solve ( [4.2| ) and obtain the expression for 77 (°) in terms of a(t), b(t) and c(t) 
or equivalently A(t),B(t) and C(t) 

rf® = -3d' 2 f Sin(i?( ^- r)) P c KM 2 dr (4.10) 
j o B 

Recall that the dependence on the internal mode amplitudes is through M, defined 
in (|3.7|). Substitution of ( f4.10| ) into equations (|4.6| -^.8|) yields a closed system for the 
internal mode amplitudes through order e 3 . The key terms in this equation come from 
certain resonances and our goal now is to show how they arise. 

The first column of terms in F involves the interactions among discrete modes. They 
do not generate frequencies contributing to any resonant forcing. The second column 
contains terms which couple the discrete bound state part to the continuum radiation 
modes. In order to identify all the resonances one has to explicitly expand out the F i2 
terms and to use equation fl4.10| ). We do not carry out the detailed computations in 
all detail here. Rather, we illustrate the key ideas and methodology by considering 
prototypical terms. A complete and rigorous implementation of these ideas in another 



nonlinear wave context is presented in |k| . 

We focus on the 0(e 2 ) term in F12 and some of its contributions to the equa- 
tion for A(t). In particular we shall consider all resonant contributions and a sample 
nonresonant contribution. By equation ( f4.6|) we must then consider the expression: 

3ie 2 d- 2 cu- 1 e- iu '<' t (Kg, Mr/ (0) ). (4.11) 

Consider r)(°'. 



r,® = -3d' 2 f Sin(jB(t T)) P c KM 2 dr = -—d~ 2 t e iB ^ P c KM 2 dr + ... 
Jo B 2iB Jo 

= ^-d~ 2 e wt f e~ Wr P c Kg 2 a 2 (r)dr + ... 
2B Jo 

= 3L e iBt [\-iBr (A 2 (T)e 2i "e T + 2\A(T)\ 2 + A 2 (T)e- 2i "* T ) P c Kg 2 dr + ... 
2B Jo ^ ' 

= ^LjBt f A 2, y-i{B-2u g )r p Kg 2 dT + ^L e iBt [* -^+2^2, ) p K tf dr + ... 

2B Jo U y 2B Jo K J 

= 4°Js + r}^ + ... (4.12) 



We have only kept two terms in the above calculation: the only term leading to a 
resonant contribution in the A equation and one (of many) nonresonant terms. 

We wish to expand r/^ using the following formula, which follows by straightforward 
integration by parts: 

e iBt /V i(B ^ )r P c a(r)dr 
Jo 

= ie iBt (B - C)~ 1 P c e iC 'a(r) - ie iBt {B - C) _1 P c a(0) 

-ie iBt C e-^ B -° T (B - Q- l P c d T a(r)dr (4.13) 
Jo 
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with a(r) = A 2 (t) P c Kg 2 , for example, and that r-derivatives of A are of order e. For 
the term i]!^) s , £ = 2u g , which for d in the range [0.54, 0.6364] lies in the continuous 
spectrum; see Table 2 and Figure 6 below. Such resonant terms are of paramount 
interest and govern energy transfer. To treat these resonant terms, an appropriate 
modification of ( |4.13| ) is required. We use the following: 
Let k = sgn(t), which is equal to +1 if t > and —1 for t < 0. 

e iBt [ t e -i{B-C)Tp ea , T ) dT = 

Jo 

ie iB \B - C + mO)- 1 P c e lC 'ct(r) - ie iBt {B -Q + m)~ 1 P c a(0) 

-ie lBt C e- i(B " c)r (5 - C + kO)" 1 ^ d T a{r)dT, (4.14) 
Jo 

where 

(B - C + kO)" 1 = ]im(B - C + z«e) -1 . (4.15) 

e|0 



Remark 4.1 (%) T/ie sense in which the operator expansions ( \4-13j ) and (U-14 ) are 
correct is in a distributional sense. That is, equality holds when multiplying both sides 
by a smooth function with spatial support which is compact and integrating both sides 
over all space. 

(2) Formula (\4-14) ^ s proved by first writing the integral on the left hand side as 



/ exp(— i(B — C + iKe)T) P c a(r) dr. 
Jo 



(4.16) 



For any e > ( {.!$ ) can be used with ( replaced by ( — it. We then pass to the limit 
as e 10. 

(3) The choice of regularization, +i0 (k = 1) for t > and —iO [k = —1) for t < is 
connected with the condition that the latter two terms in ( fj . 14\ ) consist of outgoing ra- 
diation at spatial infinity, and is therefore time-decaying in an appropriate local energy 
sense J53|/ . In this article, we take into account only the first term in ( 4-14) - Reference 
ffitf contains a fully detailed and rigorous treatment in a related context. Henceforth, 
we shall for simplicity assume t > and therefore work with the +i0 regularization. 

(4) If C does not lie in the continuous spectrum the formula j \4-H\ ) reduces to ( jJ^TQ ). 



n (o). 

Ires' 



We now continue the expansion of ?/ ) using ( |4.13| ) to study rj^J and ( |4.14| ) to study 



„(0) 
Ires 

Inr 



-d- 2 B-\B - 2u g + iO^Ait) 2 P c Kg 2 + 
\d- 2 B-\B + 2u g )- 1 A(t) 2 P c Kg 2 + ... 



Substitution of r/^ = 77^ + r}$ + ... into ( |4.11| ) we have: 



3ie 2 d- 2 uZ 1 e~ iuJst (Kg, Mrf Q) ) 
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3ie 2 d~ 2 uj g - 1 e~ iuj! ' t [Ae^ + Ae'^) (Kg 2 ,r] {0) ) + ... 
9 
2 



~ie 2 dr A u- l \A\ 2 A(Kg 2 , B~ l (B - 2oo g + iO)' 1 P c Kg 2 + ...) 



-liePd-^u^Ah-^iKff*, B~\B + 2u g y 1 P c Kg 2 ) + ... 
2 y 

= e 2 (-T 2Wg + i\ 2u)g ) \A(t)\ 2 A(t) + p(t)A 3 (t) (4.17) 

To calculate A and T, we apply a generalization to self-adjoint operators of the 
well-known distributional identity: 

lim(£ ± ie)' 1 = P.V. £~ x =F iir6{Z), (4.18) 

where 5(£) is the Dirac delta mass at £ = and P.V. denotes the principal value 
integral. 

Therefore, using (|4.18 ) we obtain: 



e 2 T 2u)g ee e 2 Y{2u g ) = ^- 2 e 2 d- i {Kg 2 , 5(B - 2cu g )Kg 2 ), (4.19) 
e 2 A 2lVg = e 2 A(2uj g ) = ~s 2 d-*(Kg 2 ,P.V.(B - 2uo g )P c Kg 2 }. (4.20) 

Remark 4.2 (1) Had we done this calculation for the B (respectively, C) equation, 
we'd have obtained as a coefficient of\B\ 2 B (respectively, \C\ 2 C ) the quantity e 2 (— T 2uJs + 
i^ 2uJs ) (respectively, e 2 (-T 2ue + iA 2u J) (2) If ( G a cont (B), then is always nonneg- 
ative and, generically, is strictly positive. It is the analogue of Fermi's golden rule 
wwhich arises in the context of the theory of spontaneous emission Apart 
from a positive constant pref actor, is the square of the Fourier transform of Kg 2 
relative to the continuous spectral part of B evaluated at Q. 

The above calculations yield the following information on the structure of the equa- 
tion for A(t): 

A t = e 2 (-T 2uJg + iA 2u)g ) \A\ 2 A + .... (4.21) 

Extensive calculations, of which the above are representative, yield a system of the 
form: 

A t = e 2 ( ai \A\ 2 + a 2 \B\ 2 + a 3 \C\ 2 ) A 

e 4 (a 4 |A| 4 + a 5 \A\ 2 \B\ 2 + a 6 \B\ 4 + a 7 \B\ 2 \C\ 2 + a 8 \C\ 4 ) A 
+ eZ A (t,A,B,C;e) 

2 (a a\1 , q l d|2 , a ly»|2 



B t = e 2 {(3\A\ 2 + (3 2 \B\ 2 + (3 3 \C\ Z ) B 

e 4 (fa\A\ 4 + I3 5 \A\ 2 \B\ 2 + (3 6 \B\ 4 + f3 7 \B\ 2 \C\ 2 + f3 8 \C\ 4 ) B 
+ eE B (t,A,B,C;e) 
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C t = e 2 {y 1 \A\ 2 + l2 \B\ 2 + l3 \C\ 2 )C 

(ia\A\ 4 + 7 5 |^| 2 |£| 2 + 7 6 |5| 4 + l 7 \B\ 2 \C\ 2 + 78 |C| 4 ) C 
+ eE c {t,A,B,C;£) (4.22) 

The coefficients atj, (3j and 7^ are, in general, complex numbers. The terms EU, 
Ea and EU involve bounded and oscillatory complex exponentials in t multiplying 
monomials in A, B, C of cubic or higher degree. Also, in order to obtain the fifth 
degree terms it is necessary to construct rj through second order in e (and therefore the 
continuous spectral part of the perturbation about the kink, e 2 rj, through e A . Coupling 
to rj is neglected as is the dynamical equation for rj. The preceding calculation gives 
«i = — T g + iA g plus a further term contributed by the 0(e 2 ) entry of Fu in ( |4.9| ). 
An exhaustive tabulation of all coefficients «j,/3j,7j would be, to put it midly, a very 
lengthy exercise. Since the principal effect we seek to illuminate is that of nonlinear 
resonant coupling on the internal mode amplitudes \A\, \ B\ and \C\, we only tabulate 
those coefficients up to the order considered which may play a role. 

For this it is convenient to introduce the notation: 

G(() = P C B-\B - C + iti)- x P D (4.23) 
Table 1A: Principal <fi 4 g-mode coefficients 



Term 


Coefficient (<Xj) 


\A\ 2 A 


-*id-*w-'{Kg\G(2u g )Kg*) 


\B\ 2 A 


-9id-*u-'-{Kgs, G(oj g + uj s )Kgs) 


\C\ Z A 


-Qid-^-^Kge, G{u g + u e )Kge) 


A^A 


-iid-^-^g^G^g*) 


\B\ A A 


-Zid-^igs*, G{u g + 2cj s )gs 2 ) 


\C\*A 


-iid- 4 uj-t(ge 2 .,G{u g + 2u e )ge 2 ) 


A\ 2 \B\ Z A 


-^id-*u- l {g"a, G(2to g + u s )g*s) 


A\' 2 \c\' 2 a 


-^d-*u>-'{g*e,G{2u> g +w e )g*e) 


\B\ 2 \C\ 2 A 


—9id~ 4 L}~ 1 (gse, G(oj g + uj s + uj e )gse) 


Table IB: Principal (j> 4 s-mode coefficients 


Term 


Coefficient (j3j) 


\B\ Z B 


-%id-*u- L (Ks A , G{2lo s )Ks z ) 


A\ 2 b 


-Md-^oj-^Kgs, G{uj g + w s )Kgs) 


\C\ Z B 


-md^oj-^Kse, G(ui s + uJ e )Kse) 


\B\ 4 B 


-|id- 4 cj s 7 1 (s M ,G(3w s )s a ) 


A^b 


-lid-^{sg\G{2u g + w s )sg') 


\C\ 4 B 


-^id-^J l {se\ G(lu s + 2uj e )se 2 ) 


A\''\B\ 2 B 


-hd-*co- 1 (gs*,G(u; g + 2u s )s 2 g) 


A\' 2 \B\' Z B 


-Ud-^-'(g S z 1 G(2cu s - u g )s 2 g) 


\B\ Z \C\' 2 B 


-\id-*u)J L {ii i e, G(2lj s + oj e )s 2 e) 


A\' z \cfB 


—9id~ 4 uj~ 1 (gse, G(uj s + uj e — u> g )gse) 


AV\c\ 2 b 


—9idr^U!~ x {gse, G{uj g + uj s + LO e )gse) 



Table 1C: Principal <f> 4 e-mode coefficients 
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Term 


Coefficient (7^) 


\C\ l C 


Ti ^ T — '. ?T ^ Ti 

~iid- i uj- L (Ke 2 ,G(2uj e )Ke 2 ) 


\B\ 2 C 


-9id- 4 'uj- 1 (Kse, G(u> s + uj e )Kse) 


\A\ 2 C 


-9id-^u}~ l {Kge, G(u g + uj e )Kge) 


|C| 4 C 


-|id- 4 w- 1 (e 3 ,G(3o; e )e 3 ) 


\B 4 G 


-^id^uj- 1 (s z e, G(2w s + w e )s 2 e) 


|^| 4 C 


-f irf-'w-^e, G(2co 5 + w e )0 a e> 


\C\ 2 \B\ 2 C 


-^rf- 4 ^- 1 (se 2 , G(w s + 2uj e )se 2 ) 


\C\ 2 \B\ 2 C 


— 7;id~ 4: uj~ L (se 2 , G(2u> e — uj s )se 2 ) 


\C\ 2 \A\ 2 C 


-Ud-^-^ge 2 , G{u g + 2u e )ge 2 ) 


\C\ 2 \A\ 2 C 


-%id-^- l (ge\ G(2uj e - oj g )ge 2 ) 


\B\ 2 \A\ 2 C 


— 4 dj~^ (gse, G(u> g + u> s + u> e )gse) 


\B\ 2 \A\ 2 C 


—9id~ 4 uj~ 1 (gse, G(uj s + uj e — u> g )gse) 



Before discussing the information contained in these tables, we note that the system 
Q4.22Q can be further simplified. Using a near-identity transformation: 

(A, B, C) = {A, B,C)+ O (\A\ 2 + \B\ 2 + \C\ 2 ) (4.24) 

the system fl42p can be transformed into a new system for A, B and C of a very 
similar form, but with the following modifications: 

• The real parts of the coefficients are the same but the imaginary parts may be 
modified. 

• The £ 2 S terms are now replaced by terms of order e 6 

We refer to the system governing A, B and C, obtained in the manner, after neglecting 
the 5 terms, dispersive normal form. 

While complicated in its details, there is a simple way to think about this normal 
form. We introduce the internal mode powers |j: 

p = \A\ 2 , Q = \B\ 2 , R = |Cf. (4.25) 

The equations for the powers are: 

P t = 2e 2 (a[P + a r 2 Q + a r 3 R) P 

+2e 4 (a r A P 2 + a\PQ + a r 6 Q 2 + a r 7 QR + a r s R 2 ) P (4.26) 
Qt = 2e 2 (P[P + (3 T 2 Q + %R) Q 

+2e 4 (f3 r 4 P 2 + (31PQ + (31Q 2 + (3 r 7 QR + (3 r s R 2 ) Q (4.27) 

Rt = 2e 2 { 1 [P + Y 2 Q + llR)R 

+2e 4 (7IP 2 + ilPQ + 1IQ 2 + Y 7 QR + 1IR 2 ) R, (4.28) 

where aj, /3J and 7J denote the real parts of the coefficients a£, /3J and 7J appearing 
in Q4.22 ). The real parts do not change under the near identity change of variables: 



7 Strictly speaking, by (4.24) P, Q and R are only approximately equal, respectively, to the Gold- 
stone, shape and edge mode powers. 
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A A, B i— > B, C \— > C, and are therefore given by the real parts of the coefficients 
displayed in Table 1. In general, as the reader may note, a contribution to the above 
mentioned internal mode power equations of the form |^4| 2m i|_g| 2rn 2|(7| 2m3 comes from a 
frequency combination miUJ g +m20J s +m^uj e landing in the band of continuous spectrum 
through a term 

(M(K, g, s, e), G{m x u g + m 2 u s + m 3 u e )M(K, g, s, e)) (4.29) 

where M is the appropriate monomial combination of the relevant spatial parts. 

As the discreteness parameter, d, varies the spectrum of the kink (the continuous 
spectrum and the number and location of the internal modes) changes. As indicated 
in Figure 6, for different values of d various integer linear combinations of the internal 
mode frequencies, the simplest of which are those appearing as arguments of G(-) in 
Table 1, may lie in the continuous spectrum. 

How does this influence the character of the normal form? 

Let ( denote one such linear combination of frequencies. By ( f4.18| ): 



G(C)=PcB-\B-C)- 1 P c , Ci<Te*nt{B) 

G(C) = P c P.V. B-\B - C)- 1 P c - i^P c 5(B - C) P c , C e a cont (B). (4.30) 

Each of the coefficients of the system listed in Table 1 is a positive multiple of an 
expression of the form —i(f, G(()f), where / is spatially localized. Therefore, if £ does 
not lie in the continuous spectrum of B, the associated coefficient of the system for the 
powers P, Q, R, a r = 9Ra, (3 r = $1/3 or Y = 9fry, will be zero. On the other hand, if £ 
does lie in the continuous spectrum of B, this coefficient will be of the form: 

Tc = -^(PJ,5(B-C) PJ) 

= -^b[/](C)| 2 , 

where Tb denotes the Fourier transform with respect to the continuous spectral part of 
the operator B. Generically, one has is strictly negative. Therefore, such resonances 
are associated with nonlinear damping of energy in the internal modes. This does not 
contradict the Hamiltonian character of the equations of motion. The A, B, C system is 
coupled to the dispersive system governing r\\ damping of the discrete mode amplitudes 
implies a transfer of energy from the discrete to continuum modes. The information 
contained in Figure 6 enables us to determine, for each d, which combinations of 
harmonics appear in the phonon band. Then Tables 1A, IB and 1C, together with 
( |4.26| - |4T28| ) give us the precise form of the internal mode power equations, from which 



we can ascertain the detailed behavior of solutions in a neighborhood of the static kink. 

In the Table 2 (below) we present the form of the internal mode power equations 
for different ranges of the parameter, d. There are numerous changes in the form of 
these equations as d varies, so for clarity, we indicate only the key transitions. These 
key transitions occur across values of d where there is a topological change in the phase 
portrait of the system governing the internal mode powers. Such changes are found 
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Figure 6: Integer linear combinations of internal mode frequencies: In all 6 panels, 
the band is indicated by dash-dotted lines. Top left shows the goldstone frequency as a 
function of d (solid) and its second (dashed) and third (dotted) harmonics for 4 . The 
top right panel follows the same sequence of symbols for the shape mode harmonics 
while the left panel of the second row does so for the edge mode. The right panel of the 
second row shows the combinations of cu g , oj s that can appear in the band. In particular: 
u g + u s :o, 2u g + u s :x, 2u s + u) g : +, 2u s — u g :*. Similarly in the bottom left panel for 
the g, e modes. In particular: uo g + u e : o, 2u g + oo e : x, 2u e + oo g : +, 2uJi — u g : *. And 
in the case of bottom right panel s, e modes are shown but so are combinations of all 3 
frequencies. In particular: u s +u e : o, 2uj s +uj e : x,2u e +u s : +,2wj — uj s : *, u s +u g +u e : 
diamonds, uj s + u) e — u g :down triangles. 
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due to a change in the nature of the set of equilibria, e.g. going from a system with 
one line of equilibria to one where there are two lines of equilibria, or due to a change 
in the number of internal modes (jump in the dimensionality of the phase portrait). 
The latter transition occurs at d — d e , d e ~ 0.82, when a point eigenvalue emerges 
from the edge of the continuous spectrum and appears as a third internal mode, e, an 
edge mode, with corresponding frequency uj e . 

Table 2: 4> 4 internal mode power equations 



Regimes of d 


Resonances 


System Form 


1 : d < 0.5398 


{2UJ S - UJg,- - ■} 


Pt = 0, Q t = -e^{PQ z } 


II : 0.5398 < d < 0.6145 

2uj g <E a cont 


2i0g 1 UJg + UJg 
2UJg UJg 


P t = -e z P{P,Q}, 

Q t = -e 2 PQ{l,e 2 {Q}} 


III : 0.6145 < d < 0.6364 


2UJg, 2UJ S - UJg 


Pt = s {P }, 
Qt = -e^PQ 2 } 


IV : 0.6364 < d < 0.6679 

2uJ g ^ <J C ont 


UJg + UJ S 
2LO s - UJg 


Pt = -e 2 {PQ}, 

Q t = -e 2 {PQ{l,e 2 {Q}}} 


V : 0.6679 < d < d e 
d e ~ 0.82 


{iUJg,AUJg,2uJ S - UJg} 
{LOg +UJ S , 2UJg +UJ S } 


Pt = -e 4 P a {l, e 2 P} - e*PQ{l, e z P}, 
Qt = s 2 PQ{l,e 2 {P,Q}} 


VI: d e <d< 0.9229 
uj e appears 


{4UJg,5UJg,2uJ S - UJg} 

ujg + uj s ,2ujg + uj s ,2uj e - uj s 

UJg + UJ e , 2UJg + U e ,UJ e + UJ S - UJg 


Pt = -e 2 P{Q{l,e 2 P},R{l,s 2 P} 1 s i P ]i } + ... 

Qt = -e 2 PQ{l,e 2 {P,Q,R}} + ... 

Rt = e 2 PR{l,s 2 {P,Q}} - e 4 QR 2 + ... 


VII : 0.9229 < d < d* 
d* - 1.2234 

2l) s 6 ffcont 


{nw 9 }5<„<38, 2uj s 

UJg +UJ S ,UJ g +UJ e ,2(jJ g +UJ S , 

2uj g + uj e , 2uj g + 2uo Sl uj g + 2uj s , 
ujg+uj s +uj e ,uj s + uj e , 2uj e — uj s 

2UJ e - UJg, UJ e +U) S - UJg, 2LO s - UJg 


P t = _^«-2{p«} + ... 

Qt = -{e 2 Q 2 } + ... 
Rt =-e 2 R{Q,P} + ... 


VIII : d*<d 

2iO e € Ucont 


nuig (5 < n < N), 2uj s ,2uj e 


Pt = -e 2n - 2 {P n } + ... 
Qt = -{s 2 Q 2 } + ... 
Rt =-{e 2 R 2 } + ... 



Remark on notation: We illustrate the notation of the table with the following 
example. Consider the regime V. This regime can be broken into three subregimes 
illustrated in the following table: 



Table 3: internal mode power equations for subregime V 



Subregime of d 


Resonances 


System Form 


a) 0.6679 < d < 0.743 


3UJg,UJg +UJ S 
2UJ S - UJg 


P t = -e 2 dPQ - e 4 C 2 P 3 , 
Qt = s 2 C 3 PQ - s 4 C 4 PQ 2 


b) 0.743 < d < 0.7622 


3Wg,4Wg 

UJg + UJ S , 2UJ S - UJg 


P t = -e 2 dPQ - e 4 C 2 P 3 - £ e C 5 P 4 , 
Q t = -e 2 C 3 PQ - e 4 C 4 PQ 2 


c) 0.7622 < d < 0.7687 


3UJg,UJg +UJ S ,2UJ S - UJg 
AUJg, 2uJg + UJ S 


Pt = -e 2 dPQ - e 4 C 2 P' d - e 4 C 6 P 2 Q - e b C 5 P 4 
Qt = -e 2 C 3 PQ - e 4 C 4 PQ 2 - s 4 C 7 P 2 Q 


c) 0.7687 <d<d e 
d e - 0.82 


UJg +UJ S ,2UJg +UJ S 
4:LOg, 2UJ S - UJg 


P t = -e 2 C x PQ - s*C 6 P 2 Q - e b C 5 P 4 
Q t = -e 2 C 3 PQ - e 4 C 4 PQ 2 - e 4 C 7 P 2 Q 



Although the details of the system form change between regimes, the topological char- 
acter of the phase portrait and therefore the qualitative nature of the solutions does 
not change; each phase portrait has P = as a stable line of equilibria. For regime 
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Ill, Table 2 is read as follows: in this regime some or all resonances occur from among 
each of the indicated sets: {Su g , 4u; 9 }, and {u g + cu s , 2u s — u g , 2u g + u s }, giving rise to 
terms in the following manner: 

ku g + lu s + mv e -> e 2n ' 2 P k Q l R m , k,l,me% + , k + I + m = n. (4.31) 

Thus the sets in curly brackets in Table 2 are to be viewed as columns of a "menu" from 
which one (or the dynamical system) chooses all or some items (resonant combinations) 
depending on the subregime of d. For a given subregime, this choice of subset gives rise 
to linear combination with nonnegative coefficients, Cj, of the corresponding monomials 
in P, Q and R in the power equations. Therefore P t ,Qt,Rt < 0. Generically , we have 
Cj > 0. 

It is straightforward to analyze the sets of equilibria and their stability for each 
of the systems in Table 2. These take the form of constant vectors with at most one 
nonzero component. These states are dynamically stable. If we take the view that the 
system of power equations determines the nonlinear dynamics near the static kink we 
anticipate that: 

• the zero solution of the internal mode power equations corresponds to the ground 
state static kink solution discrete nonlinear equation. We denote the ground state 
kink by K gs . 

• a nonzero equilibrium state with Q ^ corresponds to a time periodic solution: 

Ui ~ K g8ti + cos(u s t)si. (4.32) 

We call such a periodic solution of the full nonlinear dynamical system, which 
would have the same symmetry as the kink, a wobbling kink, which we designate 
by sW or simply W . 

• a nonzero equilibrium state with P ^ power equations corresponds to a time 
periodic solution: 

Ui ~ K gSji + cos(uj g t)gi. (4.33) 
We call such a periodic solution a g-wobbling kink. We denote this state by gW. 

• a nonzero equilibrium state with R ^ power equations corresponds to a time 
periodic solution: 

Ui ~ K gs , t + cos(u e t)gi. (4.34) 

Such a periodic solution of the full nonlinear dynamical system, would not have 
the same symmetry as the kink. We call such a periodic solution an e-wobbling 
kink. We denote this state by eW. 

Below, we present a table of the kinds of static and periodic states anticipated by 
the normal form / internal mode power equation analysis. 

Table 4: (p 4 normal form, equilibria and anticipated coherent structures 
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Regime of d 


Equilibria 


Coherent structures 


l:d< 0.5398 


{(P,0) 
{(0,0) 


P > 0} 

0>o} 


K gs , W,gW 


II-III: 0.5398 < d < 0.6364 


{(0,0) 


>o} 


K gs , W 


IV: 0.6364 < d < 0.6679 


{(0,0) 
{(^,0) 


>o} 

P> 0} 


K g8 , W, gW 


V: 0.6679 < d < d e 
d e ~ 0.82 


{(0,0) 


0>o} 


K gs , W 


VI: d e <d< 0.9229 


{(0,Q,0) :Q>0} 
{(0,0, P) : R > 0} 


K gs , W, eW 


VII: 0.9229 < d < 1.2234 
d„ ~ 1.2234 


{(0,0, P) : P > 0} 


K gs , eW 


VIII: d» < d 


{(0,0,0)} 


K 9S 



4.2 The normal form for discrete sine- Gordon 

The procedure for deriving the normal form for the internal mode amplitudes and the 
internal mode power equations presented in sections 3 and 4.1 can be applied to the 
discrete sine-Gordon equation as well. The implementation is actually simpler because 
for discrete SG there are only two internal modes: g and e; see section 2. Therefore, 
the decomposition of the solution is: 

Ui(t) = K t + ea(t) gi + eb(t)e t + s\(t), (4.35) 

where 

(g,r ] (t)} = (e,r ] (t)}=0 

PcV = V~ (9, V(t))g - (e, 7}{t))e = r](t) 

As in the previous subsection the amplitude equations for the slowly varying mod- 
ulation functions can be obtained: 

A t = e 2 (a x \A\ 2 + a 2 \B\ 2 ) A + e 4 a 3 \B\ 4 A + E A (A,B,t) 

B t = e 2 (p 1 \A\ 2 + p 2 \B\ 2 )A + e 4 p 3 \B\' l A + E B (A,B,t) (4.36) 

The coefficients otj,f3j can be evaluated along the lines detailed in the pervious section 
and are tabulated below. 



Table 5A : Principal SG g-mode coefficients 



\A\' 2 A 


-^-V,- 1 (sin Kg 2 , G{2uj g ) sin Kg 2 ) 


\A\ 4 A 


-ifW^H 1 ( cos K 9 4 , G(3uj g ) cos Kg*} 


\B\ Z A 


— jd^^ojy 1 {sixi K ge, G(cj g + lu c ) sinKge) 


|P| 4 A 


— j^yjd~ 4 uj g ~ 1 (cos K ge 2 , G(u) g + 2ui e ) cosKge 2 ) 


\A\ 2 \B\ 2 A 


— d~ 4 lu g 1 (cos Kg 2 e, G(2uj g + uu e ) cos Kg 2 e > 
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Table 5B : Discrete SG e-mode coefficients 



\B\ 2 B 


-^d-^uj- 1 (sin Ke 2 ,G(u! e ) sin Ke 2 ) 


\B\*B 


i( ^, r d-' l uj- L (coaKe- i , G(3w e ) cosKe*) 


\A\ 2 B 


— 4\d~ 4 uj~ 1 (sin K ge, G{uj g + uj e ) sin Kge) 


\A\ 4 B 


~4(3') 2 d^ 4 ^^ 1 (cos Keg 2 , G(uj e + 2uj g ) cos Keg 2 ) 


\A\ 2 \B\ 2 B 


— 4( ^yj, d^^Lu^ 1 (cos Ke 2 g, G(uj g + 2uj e ) cos Ke 2 g) 


\A\ 2 \B\' 2 B 


— jngyi d^ 4 ^ 1 (cos Ke 2 g, G(2uj e — uj g ) cos Ke 2 g) 



As with discrete 4 , the details of the internal mode power equations change as 
d varies due to the different types of resonances with the continuous spectrum which 
may occur. For discrete SG there are essentially only two regimes, and one value of the 
parameter d across which there is a topological change in the phase portraits. Since 
the detailed picture is simpler we tabulate it in greater detail. 

Figure 7 displays the variation of the internal mode frequencies (u g ,w e ), certain 

multiples of them and certain other integer linear combinations of them. As with 

discrete 4 , transitions in the structure of the normal form occur across values of d for 

which there is a change in the set of integer linear combinations which lie in the band 

of continuous spectrum. The precise normal form and power equations can be worked 

out using the coefficient Table 5 and the expression for (7(C), ( |4.30| ). Notice that only 

the leading order terms are given in these internal mode power equations. 

Table 6: SG internal mode power equations 



Regime of d 


Resonances 


System Form 


I: d < d e ~ 0.515 


None 


Pt = 


II: d e <d< .565 


2LO e - UJg, {iUJg - LO e } 


Pt = -e b {P a Q}, Qt = ~e 4 {PQ 2 } 


III: 0.565 < d < 0.65 


2ujg, 2iJ e - UJg 


Pt = -e 2 {P 2 }, Qt = -e 4 {PQ 2 } 


0.65 <d< 0.7 


2UJg,UJg + LU e 


Pt = -s 2 P{P,Q}, 
Qt = -e 2 PQ{l,e 2 Q} 


0.7 < d < 0.76 


2u g ,3ujg,ujg + uj e , 

2iO e — LUg 


Pt = -e 2 P{P,Q,e 4 P 2 }, 
Qt = -e 2 PQ{l,e 2 Q} 


0.76 < d < 0.785 


3uj g ,uj g + u e , 

2U) e — UJg 


Pt = -e 2 P{Q} - e 2 {P 2 }, 
Qt = -s 2 PQ{l,e 2 Q} 


0.785 < d < 0.8 


3UJg,4UJg 

UJg + UJ £ , 2UJ e - UJg 


Pt = -e 2 P{Q,e 2 P 2 ,e' 6 P 4 }, 
Qt = -e 2 PQ{l,e 2 PQ} 


0.8 < d < 0.847 


3UJg 7 4:Ug,2uJ e - UJg 
UJg + UJ £ , 2UJg + UJ & 


Pt = -e 2 P{Q,e 2 P 2 ,e 2 PQ,e 4 P' 6 }, 
Qt = -e 2 PQ{l,e 2 Q,e 2 P} 


0.847 <d<d* 
- 0.86 


3ujg, Aujg, hujg, 2u e — UJg 

UJg + UJ e , 2uj g + UJ e 


Pt = ~e 2 P{Q,e 2 P 2 ,£ 4 P' i ,e 2 PQ,e l3 P 4 }, 
Qt = -e 2 PQ{l,e 2 Q,e 2 P} 


IV: d*<d< 0.9 


3uj g ,4ujg,5uj gi 2uje 

UJg + UJ e , 2uj g + UJ ei 2uj e ~ UJg 


Pt = -e 2 P{Q,e 2 P 2 ,e 4 P i ,e 2 PQ}, 
Qt = s 2 Q{P,Q,£ 2 PQ,£ 2 P 2 } 


0.9 < d < 0.99 


^UJg, 5UJg, 2uJ e 

UJg + UJ e , 2uj g + UJ e , 2uj e ~ UJg 


Pt = -e 2 P{Q, e 4 P A , e 2 PQ, e 4 P 4 }, 
Qt = -e 2 Q{P,Q,s 2 PQ,e 2 P 2 } 


0.99 < d < 1.003 


5uj g , 2uj e , 2uj e - UJg 
UJg + UJ e , 2uj g + UJ e 


Pt = -e 2 P{Q,e 2 PQ,e b P 4 }, 
Qt = s 2 Q{Q,e 2 PQ,s 2 P 2 } 


1.003 < d 


nuig [n > 5), 2uj e , uj g + uj e 

2uJ g + UJ e , UJg + 2uJ e 


Pt = -e 2 P{Q,e 2 PQ,e 2n ~ 2 P n - L ,e 2 Q 2 ,...}, 
Qt = -e 2 Q{P,Q,e 4 P 2 ,e 4 PQ,...} 



The inferred coherent structures are displayed in the following table. 



Table 7: SG normal form, equilibria and anticipated coherent structures 
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Arithmetic of SG kink frequencies 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 



3 - 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

d 



Figure 7: Top panel shows the Goldstone harmonics (solid: first, dashed: second, 
thick:third and again solid: fourth). Higher harmonics of g will also be relevant, 
but are not shown here. Middle panel shows the shape mode (solid line) that lives 
close to the band edge (maximum bifurcation is of the order of 0.01) and its second 
harmonic (dashed). Finally, lower panel shows mixed combinations of harmonics (solid: 
to uj g + uj S) dashed: 2u g + u e and dotted to u g + 2u e ). The circles denote 2uj e — uj g while 
the down triangles ?>u g — u e -see note below-. In all cases the phonon band edges are 
shown by dash-dot lines. 



Regime of d 


Equilibria 


Coherent structures 


I: d<d e 


P > 


K gs , gW 


II: d e <d< 0.565 


{(P,0),(0,Q) : P > 0, Q > 0} 


K gs , W, gW 


III: 0.565 <d<d* 
d* ~ 0.86 


{(0,Q):Q>0} 


K gs , W 


IV: d* < d 


{(0,0)} 


Kgs 



5 Time periodic solutions 

5.1 Existence of wobbling kinks 

From the discussion of the previous section, we expect that the dispersive normal 
form we have derived for discrete nonlinear wave equations (in particular, discrete 
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SG and 4 ) models the dynamics in a neighborhood of the static kink solution. This 
normal form captures the effects of nonlinear resonant interactions of the internal and 
continuum mode fluctuations about the kink. Tables 4 and 7 indicate for discrete 4 
discrete SG the anticipated existence of periodic solutions which bifurcate from the 
static kink solution. In this section we prove, under suitable nonresonance hypotheses, 
the existence of such time periodic solutions. 

In tables 3 and 5 we see that these normal forms fall into two categories (i) those for 
which the origin is the only equilibrium (d sufficiently large, d > d*) (ii) d is such that 
the normal form admits one or two lines of equilibria. In case (i) we see that the zero 
solution is asymptotically stable. By ( ]3.5|) and ( |4.35| ) this suggests that the ground 
state kink is a stable attractor. At d* there is a bifurcation in the phase portrait; for 
d < d*, zero is no longer an isolated equilibrium. Rather, the normal form now has a 
line of stable equilibria passing through the origin. This then suggests an approximate 
solution of the nonlinear wave equation of the form: 



Kj + eS e? cos(m)^nj, (5.1) 



where £q denotes an internal mode Q its corresponding frequency; Bxn = &Xn an d 
e is small. If such a family of periodic solutions exists, we say this family bifurcates 
from the kink solution in the direction of the internal mode, £q. Such directions for 
bifurcation are available for d < d*. 
For d < d* are there bifurcating time periodic solutions? 

In view of the analysis of the regime d > d*, one must be cautious about reaching a 
conclusion about the large time behavior on the basis of the above approximation. In 
fact we learn from the regime d > d* that radiation damping, due to resonant coupling 
of oscillations to the continuum, may lead to the slow decay of a solution which appears 
to be time-periodic on shorter time scales. 

On the other hand, given the approximate periodic solution (|5.1|) it is certainly 
natural to attempt a perturbative (in e) construction of a true time-periodic solution. 
The approach we use is the Poincare continuation, which is standard in the persis- 



tence theory of periodic solutions of systems of ordinary differential equations ||14|| . As 
in typical perturbation expansions, one expects a hierarchy of linear inhomogeneous 
problems governing contributions to the expansion at various orders in e. Solvability 
at each order requires that one can arrange for the solution to each inhomogeneous 
problem to be time periodic of some common period. That this can be achieved in 
the problem at hand follows from the following nonresonance condition, verified in our 
numerical investigation of the spectrum of the kink: 

(NR) Assume that no integer multiple of the internal mode frequency Q lies in the 
continuous spectrum of B. 

A precise statement of the result is as follows: 

Theorem 5.1 Let K = {i£i}j e z denote a ground state kink. Assume that the linear 
operator, B, acting on the space 1 2 {TL), has a simple eigenvalue, Q, with corresponding 
eig en function, £q satisfying B£q = Qt;n- Also assume the nonresonance condition 
(NR). Then, in a neighborhood of K, there is a curve of solutions e i— > y(t; e) passing 
through the ground state kink with the following properties: 
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• There is a number eo > such that y(t; e) is defined for all e < 6q. 

• v(t;0) = K 

• y(t; e) has period 2n/u)(t) in t, where u> 2 (e) = Q 2 (l + 0(e)) is a smooth function 
with w(0) = Q. 

• y{t-e)el 2 {7L) 

• yiit-e) = Ki + ecos{uj{e)t)^ l + 0{e 2 ). 

The next result, a simple consequence of the proof of Theorem 5.1, yields the class 
of solutions which we call wobbling kinks. 

Corollary 5.1 For discrete SG, the periodic solution bifurcating from the kink in the 
direction of the spatially odd edge mode (Be = u e e) has the same symmetry as the SG 
kink around its center. For discrete 4 the periodic solution bifurcating from the kink 
in the direction of the spatially odd shape mode (Bs = u s s) has the same symmetry as 
the 4 kink. 

We seek a solution in the form: 

m = Ki + b(t)£n,i + Vi, «6 TL- (5.2) 

Substitution of ( |5.2|) into (|2.28|) and then projection of the resulting equation separately 
onto the eigenvector £q and its orthogonal complement yields the coupled system for 
the shape mode amplitude and the radiation components 

d 2 b + n 2 b = -d-'i^VUn + ^n + r]) 2 ) (5.4) 
d 2 V + B 2 ri = -d-'P.Vlb^ + ^n + r]) 2 , (5.5) 



where 



V[M] = [ 1 (l-9)V"'(K + 6M)d9. (5.6) 
Jo 



Remark 5.1 It is possible to formulate the proof of Theorem 5.1, as essentially a 
consequence results in (Chapter 14, Theorem 2.1), applied to (\5.Jj^).d[ ); a simple 
generalization of the result to an infinite dimensional setting is required. In this context, 
hypothesis (NR), corresponds to the hypothesis that 1 be a simple Floquet multiplier 
of the unperturbed linear problem obtained by setting the right hand sides of ^5.^\5.5\ ) 
equal to zero. In this paper we present a direct proof of Theorem 5.1 in order to make 
our study elementary and self-contained. 

8 Here we employ the Taylor expansion for the function W = V: 

W(K + M)=W(K)+W'(K)M+ [ (1 - 6)W"(K + 6M)d6M 2 (5.3) 

Jo 
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We shall be considering small perturbations of the static kink, and therefore introduce 
rescaled internal mode and radiation components: 

b = eb 1 , r} = er} 1 , (5.7) 

where e is a small parameter. Since the nonlinear wave equations considered are au- 
tonomous, we can't a priori specify the period, we leave the period unspecified at 
this stage and introduce a new time variable, r, with respect to which the sought-for 
periodic solution is 2n periodic. Let 

t = uj(e)~ l T, co(0) = tt, (5.8) 

and 

&i(*) = /Mr), /3i(r + 27r)=A(T). (5.9) 
Then, equations (|5.4H5.5| ) become: 

(u 2 (e)d 2 + Q 2 )f3 1 (r) = e^Cft.ifrje) 

W^e) = -d- 2 (^,V[e(P^ n + r ]l )}(p 1 ^ + r ]l ) 2 ) 

(5.10) 

(io 2 (e)d 2 + B 2 ) m (T,t) = dFM^e) 

= -d- 2 P c V[e(^n + ViWiCn + m) 2 (5.11) 

Let 

u 2 (e) = n 2 (l + ea(e)) } (5.12) 
where a(e) is to be determined. Then ( |5.1(J| ) becomes 

(d 2 + l) p 1 = \ > (e)A + Q-^faTK e)]. (5.13) 
v 7 1 + eo~{e ) 

We extract the leading order behavior by defining: 

ft = cos(r) + eft, (5.14) 

where 

(d T 2 + l)ft = - ^[^(e) cos(r) + e <x(e)ft + fi" 2 ^ (cos(r) + eft, 77a; e) (5.15) 

1 + ecr(e) 

By the Fredholm alternative, equation (|5.15|) has an even 2ir periodic solution if and 
only if 

P cos (a cos(r) + eaf3 2 + ^Tp) = 0. (5.16) 

Here, 

-Pcos 9 = 7T _1 / coa(fi)g(ji)dfi (5.17) 

JO 
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is the projection onto the even 2ir periodic null space of d 2 + 1. The orthogonality 
constraint ( |5.16| ) can be rewritten as 

na + to / cos(t)/3 2 (t)g?t + uo~ 2 / cos(r)JFg (cos(r) + e/3 2 (r), rji; e) dr = (5.18) 
Jo Jo 

We have therefore reformulated the problem of finding a periodic solution of the discrete 
nonlinear wave equation, or equivalently (5.4- 5~5|) , as the problem of finding 2n periodic 
solutions in r, (/3 2 (e), 771(e), <r(e)), of the system: 



0. 



(5.19) 



where F = (F 1; F 2 , F 3 ) , is defined by: 
F^farn^e) = (d 2 T + l)p 2 



i-i 



-(1+6(7 



a cos(r) + ecr/3 2 + ^ .F/j (cos(r) + e/? 2 , ?7i; e) 



(5.20) 



ft 2 (l + ea)d 2 T + 5 2 ) t^(t, i) - e^(cos(r) + e/3 2 , r/ i; e) 



(5.21) 



2tt 



7ra + ea- / cos(r)/3 2 <ir 
o 



F3(02,77i,a;e) = -(1 + ea) 

f 27T 

/o 

We view F as a mapping of ((, e) G X \— > F((, e) G 3^- Here, A? and y are defined 



/■2tt 

+ fT 2 / cos(r)J r s (cos(r) + e/3 2 (r), rji; e)dr 
Jo 



(5.22) 



by: 



X 



y 



C = ?y, o") such that 
P = P(t) G i? 2 , even, and 27r periodic 
T) = T](t, ■) E H 2 even, and 27r periodic in r 
with values in the space of 1 2 {7L) functions and a G IR 
( = (/?, T], p) such that 
P = j3(r) G L 2 , even, and 27r periodic 
r] = r)(r, ■) G L 2 even, and 27r periodic in r 
with values in the space of l 2 (7L) functions and 



P 



2?r 



cos(/i)<i/i. 



We find a particular G A" for which F^ ^; 0) = 0, and then seek to construct 
curve of solutions e i— > C( e ); C(0) = C j f° r an sufficiently small e using the implicit 
function theorem fill . To find C we set e = and consider the system F((; 0) = 0. 



Taking 771 = n\ = and a = cr( ) = 0, we find that /3 2 satisfies the equation: 

(d 2 + l) /3f = -d- 2 ir 2 <^, V[0]eg> cos 2 (r), 



which has the solution 

M 0) 



]rtT 2 (^V[0]£ 2 > (l-±cos(2r) 



(5.23) 
(5.24) 
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Thus, C (0) = (/92 0) ,0,0) G X satisfies F(C {0) ;0) = 0. We shall now use the implicit 
function theorem to continue this solution to show that this solution deforms uniquely 
to nearby solutions for e sufficiently small and nonzero. 

To apply the implicit function theorem, it suffices to check that d^F((^;0) is 
bounded and invertible. A computation yields: 



^(C (0) ;0) 



d T 2 + l -cos(r) 

tt 2 d 2 + B 2 

-7T 



(5.25) 



Invertibility of d c F((^;0) can be shown by solving the system of inhomogeneous equa- 
tions: 



where 



The third equation implies: 





' 5(3 2 - 




' B ' 




Srji 


, R = 


E 




5a _ 




E 



d c F -5( = R 



and £ 



2?r 



cos(r)B(r)dr. 



<S<7 = -7r _1 E. 

Substitution into the first equation yields the equation an <5/3 2 : 



(fi + \)5I3 2 = B- it- 1 cos(n) B(fi)dfi cos(r), 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



which has an even 2n periodic solution. Finally the second equation can be rewritten 

as 



(Q 2 d T 2 + B^Srji = E = Y, cos(nr)e n 

n=0 



(5.30) 



which can then be solved by Fourier series: Srji = Y,T=o9nCos(riT). We obtain 

g n = (B 2 -n 2 Sl 2 )- 1 e n . (5.31) 

By (NR), e n is well-defined for each n. Note also that there is a strictly positive 
minimum distance of the set {nfl} to a(B). Therefore, 



lbn||p(z) < dist ((nfi) 2 ,o-( J B 2 )) ||-E n |h 2 (z), 
from which we have boundedness of the inverse: {Vt 2 d 2 + B 2 ) 1 : 

||^i||h 2 (R;/ 2 (Z)) < C||-E||l 2 (R;/ 2 (Z))- 



(5.32) 



(5.33) 



By the implicit function theorem there is a number e > such that for e < e the 
nonlinear wave equation has a unique periodic solution of the form: 



yi(r, e) = Ki + (e cos(r) + e 2 p 2 (r))^n,i + 071(1, r; e) 



(5.34) 



38 



with r = uj(e)t, ou 2 (e) = Q 2 (l + ecr(e)) and /32,Vi ^ n periodic functions of r. Note that 
i]i = 0(e), by the equation F 2 = 0, so that we have (|5.1| ). This completes the proof of 
Theorem 5.1. 

To prove Corollary 5.1, we note that by hypothesis, the direction of bifurcation is 
spatially odd. It is simple to check that the entire proof goes through with the spaces X 
and y additionally constrained to consist of functions r], such that rj(-, —i) = —T}(-, i). 



5.2 Do quasiperiodic solutions bifurcate from the kink? 

If we attempt using the above method of proof to construct quasiperiodic solutions, we 
do not succeed. To be specific, suppose we seek to construct a quasiperiodic solution 
which is generated by the two internal modes £m and Formally, we seek Ui(t) = 
u(i,t) in the form: 



u(i,t) = K(i) + e^aj6fy(z) cos(f^) + er]i(i,t) 
i=i 



Substitution into the nonlinear wave equation gives 



(d 2 + R 



2 



<Xj£n s cos(ftjt) + 771 



Expansion of i]i(i,t) as a power series in e, 



(5.35) 



3=0 

leads to a hierarchy of inhomogeneous equations of the form: 

(d 2 + B 2 ) v { j) =S( j \t) 



(5.36) 



(5.37) 



where as j increases (t) contains a finite sum of terms with time-frequencies of the 



form riiOi + n2^2 with |m| < , |na| < N2" and N[ J> , N^' increasing with j. Thus 
it is natural to solve for each rj^ as a truncated multiple Fourier series: 



U) 



7-00 TuCO i, 



E 



j(niQi+n 2 02)^ „(j) 



»2' 



(5.38) 



\ni\<N[ j) ,n 2 <N^ j) 



Substitution into (|5.37| ) yields: 

(b 2 - (nxOi + n 2 fi 2 ) 2 ) ^n 2 = <?£ 



n 2 ' 



(5.39) 
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where denotes the term in S^(t) which is proportional to exp(i(niQi + n 2 Vt 2 )t). 

In order to ensure the general solvability of (|5.39|) we need that for any n\,n 2 G 7L, 

n 1 Q 1 + n 2 n 2 i ±a(B). (5.40) 

Although it is nongeneric for a frequency in the point spectrum (internal mode fre- 
quency) to be hit, since the continuous spectrum (phonon band) is an interval, gener- 
ically one has n\Vt\ + n 2 Q 2 £ o'cont(B), for infinitely many choices of n x ,n 2 . For such 
choices the operator on the right hand side of ( |5.39| ) is not invertible and these resonant 
frequencies are an obstruction to solvability. A closer look at the set of resonances and 
their contribution on the dispersive normal form, would give insight into the lifetime 
of such eventually decaying quasiperiodic oscillations. 

6 Large time behavior in a neighborhood of K gs 

In this section we combine the normal form analysis of section 4, and the existence 
theory for time periodic solutions of section 5 with numerical simulation to get a more 
detailed picture of the large time behavior of discrete <p A and SG in a neighborhood 
of the ground state kink. We shall make repeated use of Tables 2 and 4 for 4 and 
of Tables 6 and 7 for SG in which the normal form / power equations and coherent 
structures are tabulated. In interpreting these tables we recall that e, introduced in 
( |3.5| ), measures the size of the component of the perturbation about the kink in the 
internal mode subspace. In full simulations of the evolution equation e is typically of 
order 10 _1 , and therefore some of the decay phenomena anticipated by our analysis, 
occur on very large time scales (e.g. r ~ e~ 2 , e~ A or longer) and are difficult to simulate 
accurately. 

In each of the various d— ranges we proceed as follows: we mention the relevant 
coherent structures (the high energy kink is omitted because it is unstable). These are 
all anticipated by the normal form analysis and their existence is established rigorously 
in section 5. We then discuss the numerically observed large time behavior for different 
classes of initial conditions. This gives evidence of the attracting nature of the periodic 
solutions constructed in section 5 in various regimes of the discreteness parameter d. 

As mentioned in the introduction, this is related to the final stages of the evolution 
of a propagating kink, pinned to a particular lattice site, and its damped oscillation to 
an asymptotic state within the Peierls-Nabarro potential. 

6.1 Discrete 4 

Regimes II, III and V: The coherent structures of interest are: K gs and W . Nu- 
merical simulations show that the wobbling kink, W, is a local attractor. 
Regimes I,IV: The coherent structures of interest are: K gs , W and gW. 

Experiment 1 (data given by K gs plus a small multiple of the Goldstone (even) 
g-mode): solution appears to approach gW. 

Experiment 2 (data given by K gs plus a small multiple of the shape (odd) s-mode): 
solution appears to approach W. 
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Figure 8: Simulation of discrete 4 with d = 0.8 in Regime IV. Upper panel: u(t) — K gs 
is seen to be approximately equal to the (spatially odd) shape mode. Lower panel: 
Projection of u(t) — K gs onto the shape mode, s; principal frequency of oscillation is 
approximately u s 
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Experiment 3 (data given by a general small perturbations of K gs ): gW 
appear to have basins of attraction. Although an exact determination of this 
not analytically tractable, its projection onto the internal mode subspace, the 
{g, s}, can be approximated using the internal mode power equations and the 
information on coefficients from Table 1. In particular, we find 



dQ 
dP 



and W 
basin is 
span of 
explicit 



(6.41) 



from which one gets the prediction that for data with 



|E(0)| 2 <^|A(0)| 2 , 



LO, 



(6.42) 



solutions asymptotically approach a g-wobbler, gW, while for 



|£(0)| 2 >^|A(0)| 2 , 



(6.43) 



solutions asymptotically approach a wobbler, W . 
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Figure 9: Simulation of discrete 4 with d = 0.65 in Regime IV. Upper panel: u(t) — 
K gs is seen to be approximately equal to the (spatially even) Goldstone mode, g. 
Lower panel: Projection of u(t) — K gs onto the Goldstone mode; principal frequency 
of oscillation is approximately u g . 

Regime VI: The phase portrait of the normal form jumps from dimension 2 to di- 
mension 3 due to the appearence of a (spatially even) edge mode. In addition to the 
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ground state kink, K gs , there are time-periodic (wobbling solutions) W and eW. Nu- 
merical simulations indicate the presence of very long lived quasiperiodic oscillations 
about the kink. These appear to be oscillations of the form represented in the first two 



terms of (|5.35|) . Strong evidence of the eventual decay of such oscillations can be seen 
as follows: 

• In section 5.2 we have shown that an attempt to construct quasiperiodic solutions 
breaks down due to a high order resonance, i.e. niu s + n 2 uj e £ a cont (B), where 
| rii| + |n 2 | is large. In this particular case, we have — u s + 2oj e £ a cont (B); see 
Table 2. 

• Such resonances correspond to obstructions in the power equations to equilibrium 
solutions of the form: (0, Q, R) with both Q and R nonzero. The obstructing 
term is the term —e A QR 2 in the R equation. 

• Linearization of the power equations about any equilibrium point (0, Q eq , 0) or 
(0, 0, R eq ) shows that each line of equilibria is asymptotically stable, correspond- 
ing to the conclusion that quasiperiodic oscillations will damp with a resulting 
asymptotically periodic state, W or eW, depending on initial conditions. The 
time scale of decay of these oscillations, r, is set by the order in e of the obstruct- 
ing terms, e.g. r ~ e~ 4 . 

Regime VII: The coherent structures of interest are: K gs and eW. 
Experiment 1: For data which is an odd perturbation of K gs . K gs is the attractor. 
Experiment 2: For general data, eW appears to be a local attractor. 
Regime VIII: The coherent structure of interest is K gs . 

Experiment 1 (data given by K gs plus an small perturbation in the direction of the 
(odd) shape mode): approach to K gs at a rate 0(t~^). 

Experiment 2 (data given by K gs plus an small perturbation in the direction of the 
(even) edge mode): approach to K gs at a rate 0{t~*). 

Experiment 3 (data given by K gs plus an small perturbation in the direction of the 
(even) Goldstone mode): perturbation from K gs appears to decay but at a different 
rate and on much longer time scales than in Experiments 1 and 2. 



6.2 Discrete sine-Gordon 

Following the model of the previous subsection, we indicate relevant coherent structures 
and briefly discuss the dynamics in a neighborhood of the ground state kink for the 
discrete sine-Gordon equation. 

Regime I: The coherent structures of interest are: K gs and gW. For general data, 
gW appears to be a local attractor. 

Regime II: The coherent structures of interest are: K gs , W and gW. The dynamical 
behavior is as, for example, in regime IV for the 4 model. 

Regime III: The coherent structures of interest are K gs and W. Numerical simulations 
show that W is a local attractor. 
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Figure 10: Simulation of discrete 4 with d — 1.1 in Regime IV. Upper panel: u(t) — 
K gs is seen to be approximately equal to the (spatially even) edge mode, e. Lower 
panel: Projection of u(t) — K gs onto the edge mode; principal frequency of oscillation 
is approximately uj e . 
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Regime IV: For all d > d*, there are positive integer multiples of both u g and u e 
which fall in the phonon band, a cont (B). Therefore, hypothesis (NR) of section 5 is 
not satisfied for either internal mode, and the existence theory of periodic solutions 
breaks down. A general perturbation about K gs will excite the internal modes and 
these resonances are responsible for transfer of energy from the internal modes to K gs 
and to radiation modes. The asymptotic state is observed to the ground state kink, 
with the details of the damped oscillatory approach to it, being predicted by the various 
d— dependent normal forms. 

Within each subregime of regime IV, the power equations are, up to terms which 
can be shown to be negligible for large t, of the form 

p t = - e 2n P n 
Q t = -e 2 PQ. 

It follows that 

P(t) = P (l + ( n - l)P n -V n tr^ (6.44) 

and that Q(t) is decaying at exponentially. These predicted power laws for P(t) (on a 
time scale of order e~ 2n , and the implied very rapid decay of Q are, to a good degree of 
approximation, observed in numerical simulations of the discrete sine-Gordon equation. 
In particular, using the form of the near-identity transformation (4.24), we obtain for 
the decay rates (of initial conditions with general data): 

• 0.86 ~ d* < d < 0.9: n = 3, P(t) = 0(H) -> |a(f)|, |6(*)| ~ 0(H) 

• 0.9 < d < 0.99: n = 4, P(t) = 0(H), -> \a(t)\, \b(t)\ ~ 0(H) 

• 0.99 < d < 1.003: n = 5, P(t) = 0(t~*), -»• |a(t)|, |6(*)| ~ O(rs) 

Note, however, that for purely odd initial perturbations around the ground state 
kink in this regime, the edge mode projection decays as t -1 / 2 , as dictated by the —e 2 Q 2 
term. 



7 Summary and open questions 

Summary: In this paper we have considered the large time behavior of solutions to 
discrete nonlinear wave equations, particularly the discrete sine-Gordon and discrete 
4 models, for initial conditions which are a small perturbation of a stable (ground 
state) kink. 

(1) We have proved in the regime where the discreteness parameter is sufficiently 
small, corresponding to sufficiently large lattice spacing or weak coupling of neighboring 
oscillators, the existence of various classes of time periodic solutions. These include the, 
in the case of the 4 the wobbling kinks (W), which have the same spatial symmetry 
as the kink, and which were anticipated in previous numerical studies as well as e- 
wobblers (eW) and g- wobblers (gW), periodic solutions which do not respect the spatial 
symmetry of the kink. In the SG case, it is the eW's that respect the spatial symmetry 
while the gWs do not. 
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(2) We have used the methods of scattering theory and ( Hamiltonian dispersive) nor- 
mal forms to study the final stages of pinning of a kink to a particular lattice site. For 
large values of the discreteness parameter ( "near" the continuum regime) , the asymp- 
totic state is a static kink. In contrast, for small values of the discreteness parameter, 
the above time periodic states can be attracting orbits. This is in sharp contrast to 
the behavior of solutions for the corresponding continuum equation in which a (pos- 
sibly moving) kink is the attractor. Our analysis makes clear the relation of broken 
(Lorentz) invariance to these contrasting dynamics. Our results give a systematic clar- 
ification of the physicist's heuristic picture of the dynamics of the center of mass of the 
kink as the effective (radiation-) damped motion of a massive particle in the Peierls- 
Nabarro potential. The approach to asymptotic state is via a slow damped (periodic 
or quasiperiodic) oscillation. The damping of the oscillation can indeed be very, very 
slow, and the lifetime is deducible from the normal form. 

The methods we use are very systematic and general, and apply to many situa- 



tions |51|, [52], [53| where the dynamical system can be viewed as the interaction of two 
subsystems: a finite dimensional subsystem, here governing the kink and its internal 
oscillations, and an infinite dimensional dynamical system, here governing dispersive 
radiation. In particular, one problem of related interest in which these methods can be 
readily generalized is the behavior of pulse-like breathing modes in the discrete non- 
linear Schrodinger equation iipi jt = —k5 2 ipi — \ipi\ 2 ipi- Using the monochromatic gauge 
symmetry of this equation, we can convert the time periodic problem of the breathing 
solutions into a static problem by looking for solutions of the form ipi = exp(\/— lu>t)ui 
and studying their stability in the frame rotating with the same frequency u (the so- 
called rotating wave approximation). The excitation and nonlinear Lyapunov stability 
of such states (see for example footnote 6 in section 2) is established in This notion 
of stability, being defined in terms of conserved integrals, is insensitive to the radiative 
behavior and the asymptotic approach to a ground state. For a recent account of the 



(numerical) construction and stability of such modes, see for instance [32]. Then, once 
the problem has been posed on the rotating wave frame, it has become a static problem 
amenable to the same techniques for the study of dispersive waves and the radiation 
losses of the pulse-like coherent structure as the kink-like structures studies in this 
work. Exponentially small effects in the solutions for this class of systems arising for 
problems with a perturbed form of nonlinearity have been explored using perturba- 
tion theory techniques f45| , f46|[ . These techniques are applicable to a restricted regime 
of parameter space because of the adiabatic approximation they entail. However, to 
capture, primarily, radiation effects caused by the discretization in the dynamics of 
continuum-like coherent structures our technique is obviously most appropriate (and 
clearly by no means restricted in its applicability) as illustrated by the analysis given 
above. What's more, the phenomena present in such an analysis are in general of the 
power-law type rather than exponential. 

Open questions: There are many related open problems of which we now mention 
several. 

• Although the normal form analysis anticipates the existence and stability proper- 
ties of various time-periodic solutions, we only have a rigorous theory concerning 
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their existence. Linear stability theory would require an infinite dimensional Flo- 
quet analysis, and nonlinear stability theory would require the analogue of our 
normal form in a neighborhood of a periodic, rather than static, solution. Numer- 
ically, the asymptotic stability of the time periodic structures can be concluded 
from their persistence for long time integrations (to the extent of our compu- 
tational power). The periodic structure construction and stability can also be 
performed using a limit-cycle type of technique similar to the one considered in 
[0. [10]. The rigorous stability theory however is still a challenging open mathe- 
matical problem. 



A rigorous infinite time theory, as in |p3[ , for discrete nonlinear wave equations 
with kinks is a challenging open question. The problem is challenging even in 
the continuum case, where the kink is the sole attractor. In this case, the normal 
form yields the correct asymptotic rate of decay |39|, f49|. However, due to the 



strength of the interaction (in space dimension one) of the dispersive and bound 
state components, the deviation of the full solution, u(t), from the kink plus its 
internal mode modulations is not asymptotically free, i.e. not a free dispersive 
solution of the linear Klein-Gordon equation. Rather, it requires a solution- 



dependent logarithmic in time corrected phase H49| , characteristic of long range 
scattering problems. 

Of interest would be more detailed long time simulations which would further 
elucidate the large time dynamics and assist in mapping out basins of attraction 
for various periodic states. 

Our study gives a detailed picture of the final stages of the pinning of a kink to a 
particular lattice site. However, the early stages of a kink plus large perturbation 
propagating in a lattice have some of the same features of the regime we con- 
sider. Namely [[48| , as the kink's center of mass moves under the influence of the 
Peierls-Nabarro potential, it executes repeated acceleration and slowing with an 
approximate PN frequency. This oscillation appears to resonate with continuous 
radiation modes, resulting in an emission of energy, deceleration of the kink and 
eventual capture of the kink by a particular lattice site. It would be of interest 
to extend and apply the ideas of the current paper to this problem. 

In this paper we have analyzed the problem of radiative effects of discreteness 
for the 2 — 7r kinks within the P-N barrier (for a recent review on the effects of 
discreteness as well as the physical applications of a number of models similar 
to the ones considered here see ||). It is well-known however that in the SG 
lattices multiple kinks (2nir kinks in general) can also exist. For these kinks first 
mentioned in |^8|] the static stability picture was analyzed in 0. There are a 
number of open problems regarding the behavior of such modes: 

1. Their motion prior to pinning is only very slightly effected by radiative phe- 
nomena |48|. This phenomenon hasn't been accounted for satisfactorily to 



the best of our knowledge to date. In fact, this forms a part of a more 
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general problem concerning the motion of coherent structures in lattices. 
It is known that under certain conditions and for potential different than 



the ones considered here |T^, |26|] there exist lattice systems with travelling 
wave solutions. Hence, it would be very interesting to elucidate the general 
conditions under which such lattice systems support travelling wave solu- 
tions (or equivalently under which conditions solutions to the advance-delay 
equations of the travelling wave frame exist). 

2. When trapped (eventually) in the PN barrier these moving 4 — it kinks also 
radiate in a way that can be captured by the analysis of this paper. The 
problem there will have more (as shown in |H) internal modes bifurcating 
from the continuum and hence will be more complicated but it can still be 
analyzed in the way presented above. 
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